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ABSTRACT 

We present low- and high-resolution Spitzer/IRS spectra, supplemented by 
IRAC and MIPS measurements, of 22 blue compact dwarf (BCD) galaxies. The 
BCD sample spans a wide range in oxygen abundance [12+log(0/H) between 7.4 
and 8.3], and hardness of the interstellar radiation field (ISRF). The IRS spectra 
provide us with a rich set of diagnostics to probe the physics of star and dust 
formation in very low-metallicity environments. We find that metal-poor BCDs 
have harder ionizing radiation than metal-rich galaxies: [0 iv] emission is ^ 4 
times as common as [Fell] emission. They also have a more intense ISRF, as 
indicated by the 71 to 160 /im luminosity ratio. Two-thirds of the sample (15 
BCDs) show PAH features, although the fraction of PAH emission normalized 
to the total infrared (IR) luminosity is considerably smaller in metal-poor BCDs 
(~0.5% ) than in metal-rich star-forming galaxies (~10%). We find several lines 
of evidence for a deficit of small PAH carriers at low metallicity, and attribute this 
to destruction by a hard, intense ISRF, only indirectly linked to metal abundance. 
Our IRS spectra reveal a variety of H 2 rotational lines, and more than a third of 
the objects in our sample (8 BCDs) have ^ 3a detections in one or more of the 
four lowest-order transitions. The warm gas masses in the BCDs range from 10 3 
to 1O 8 M , and can be comparable to the neutral hydrogen gas mass; relative to 
their total IR luminosities, some BCDs contain more H 2 than SINGS galaxies. 
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Introduction 



First IRAS and ISO, then SCUBA and COBE, and most recently Spitzer have convinc- 
ingly shown that most of the star formation in the universe is enshrouded in dust. Large 
populations of infrared-luminous galaxies at z £ 1.5 contribute to about 70-80% of the far- 
infrared and 30% of the sub-millimeter backgrounds, and are probably responsible for most 



of th e star-formation activity at high redshift (e.g., IChary fc Elbazl 12001 



Le Floc'h et al. 



20051 ; iDole et al.l 120061 ). Indeed, half the e nergy and most of the photons in the univers e 



come from the infrared spectral region (e.g.. lHauser fc Dwekll2001t iFranceschini et al.ll2008l ). 



That dust is so prominent in the high-redshift universe may appear surprising, since 
it has been assumed that dust would be absent in low-metallicity environments. However, 
recent theoretical developments challenge the assumption that dust is absent in metal-poor 
gas. Large amounts of dust can be created on short timescales by supe rnovae (SNe) and 



Asymptotic Giant Branch s tars which evolve in a metal- free ISM (e . g.. iTodini fc Ferrara 



2001 



Schneider et al. 2004; Nozawa et al. 2007; Bianchi fc Schneider 2007; Valiante et al. 



20091 ). Observationally, the copious dust emission observed at very high redshifts (z ^ 6 ) 
in quasars implies that indeed dust can form rapidly at early epochs (jBertoldi et al.ll2003l ). 
Although luminous quasars can hardly be considered as typic al examples of star-for ming; 
environments, observations of sub-millim eter galaxies (SM Gs, IChapman et al.l 120051 ) and 
dust-obscured galaxy populations (DOGs. iDey et al.l 120081 ) also suggest that early star for- 
mation episodes can also be very intense and relatively brief. However, exactly how these 
massive starbursts occur and evolve is not yet clear. The short interval in which star forma- 
tion and the ensuing chemical enrichment and dust formation convert a dust-free metal-free 
environment to a dusty metal-rich one by redshift ~6 is as yet unobserved, and studies of 
such transitions remain a major observational challenge. 

The Local Universe is home to star-forming dwarf galaxies that are of much lower metal- 
licity than those observed thus far at high redshift. Over the last thirty years or so, only a 
handful of blue compact dwarf (BCD) galaxies have been discovered with 12+log(0/H)~7.2. 
Despite intensive searche s, only one BCD, SBS 0335— 052 W, is known to host Hll regions 
with 12+log(0/H)<7.1 (llzotov et al.l 120091 ). Even at these extremely low metallicities, 
BCD spectral energy distributions (SEDs) can be dominated by infrared dust reprocess- 
ing of ultraviolet radiation from young massive stars (e.g., S BS 0335— 052 E and IZwl8: 



Thuan. Sauvage. fc Maddenlll999l ; iHouck et al.ll2004bl ; IWu et al.ll2007[ ). Because such galax- 
ies are chemically unevolved, they can provide a window on primordial galaxy formation and 
evolution. This is, in some sense, a "local" approach to a cosmological problem. If we can 
study the properties of a metal-poor ISM and its constituents locally, we may be able to 
better understand the high-redshift transition from metal-free Population III (Pop III) stars 
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to the chemically evolved massive galaxies typical of the current epoch. 



We adopt here this "local" approach. With the aim of characterizing the properties 
of star and dust formation in a metal-poor ISM, we have used the Spitzer Space Telescope 
to study a sample of nearby low-metallicity star-forming dwarf galaxies. This is the third 
in a series of papers; the first two papers h ave focused on tw o particular BCDs, Haro 3, 
the highest metallicity object in our sample ( jHunt et al.l 120061). and Mr k 996, a BCD with 
an unusually dense and compact nuclear Hll region (IThuan et al.ll2008l ). Here, we present 
for the first time the entire sample of 23 BCDs, and report on its spectral properties, as 
derived from IRS observation^. The sample has been carefully chosen to span a wide 
range in oxygen abundance: 12+log(0/H) varies between 7.4 and 8.3, with a median of 7.9. 
About a quarter of the objects in the sample are in the eXtremely Met al-Deficient (XMD) 
or Extremely Metal-Poo r Galaxy (EMPG) range (12+log(0/H) < 7.65 IPustilnik fc Martin 
20071 ; iBrown et al.ll2008l ). Thus our sample comprises the largest number of XMDs/EMPGs 
to date with good signal-to-noise IRS spectra. 

In §0 we describe the sample selection criteria, together with the Spitzer observations 
and the data reduction method s. In &j we derive fluxes of lines and features in the spectra 
by fitting them with PAHFIT (ISmith et al.l 120071 ). We present and discuss our results for 
the infrared (IR) fine-structure (FS) lines in §HJ the aromatic features in §0 and in §61 the 
molecular lines. Our conclusions are summarized in £0 



2. Sample Definition, Observations, and Data Reduction 



One of the most important factors in shaping the SED at low metallicity is the hardness 
and intensity of the interstellar radiation field (ISRF). The ISRF hardness and intensity 
govern dust properties such as the grain size distribution and temperature. Hard and intense 
ISRFs may destroy small grains all together, thus suppressing the aromatic or polycyclic 
arom atic hydro c arbon features (PAHs) which are typi cal of metal-rich starburst galaxies 



(e.g., iVoitl Il992l : Madden et al.ll2006l : ISmith et al.l 120071 ) . In the post- ISO and post- Spitzer 
eras, ratios of IR lines of n eon have often been u s ed to quantify th e hardness of the ISRF in 
star-forming galaxies (e.g., iThornley et al.ll2000t iDale et al.l 120061 ) . However, since we have 
selected the objects in our sample on the basis of their optical spectra, we have quantified 
ISRF hardness optically, with the ratio of the nebular Hell A 4686 emission line relative 
to H/3. The Hell line, with an ionization potential of 54.4 eV, is a good indicator of a 



1 We present IRS results for only 22 BCDs, since one, SBS 0940+544, was observed in a Guaranteed Time 
Program. 
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hard ultraviolet (UV) radiation field. Nebular Hell lin es are predicted to b e very strong 
in metal-free (or nearly metal- free) primordial galaxies (Schaerer 2002 , 20031). although no 



Hell-emitting high-redshift galaxy (z £ 5) has yet been discovered (e.g.. 



Nagao et all 120081 ). 



However, some low-redshift low-met allicity BCDs are known to show strong nebular 
Hell A4686 emission with H ell/H/3 £ 3% ( jGuseva. Izotov fc Thuan I I200CM ; llzotov fc Thuan 
2004J ; iThuan fc Izotovll2005l ). We have selected 23 objects so that they span a wide range in 
Hell/H/3. Roughly 1/3 of the sample consists of objects having strong Hell emission, with 
Hell/H/3 ^ 2%; another third have intermediate Hell emission with Hell/H/3~l-2%, and the 
remaining third has weak or no Hell emission, with Hell/H/3 £ 1%. This wide range of Hell 
strength in our sample also p rovides a large range in metallicity as there is a loose correlation 
between the two quantities (IThuan fc Izotovl 120051 ). The general properties of the sample 
are given in Table [U 

The observations presented in this paper are part of our GO Spitzer program (PID 
3139) for 23 BCDs. In the course of this progr am, we have obtain ed IRS spectra in the 
low- and high-resolution mo dules (SL, SH, LH iHouck et al.l l2004a[ ). I RAC images at 4. 5 
and 8/im (Fazio et al.l 120041 ) . and MIPS images at 24, 71, and 160 /im flRieke et al.l 120041 ) . 
For all instruments, the data reduction starts with the Basic Calibrated Data (bed). The 
corresponding masks (the DCE masks), furnished by the Spitzer Science Center pipeline, 
were then used to flag potential spurious features in the images, such as strong radiation 
hits, saturated pixels, or nonexistent /corrupted data. 



2.1. IRS spectroscopy 

Spectroscopy was performed in the staring mode, in both orders of the Short Low- 
resol ution module (SL1. SL2) and with the Short and Long High-resolution modules (SH, 



LH) ( IHouck et al.ll2004af ). Thus, we have obtained low-resolution spectra from 5.2 to 14.5 /im 
(R~64-128), and high-resolution spectra from 9.6 to 37.2/xm (R~600). We have used differ- 
ent acquisition schemes, depending on the brightness of the source as measured by its scaled 
H/3 flux. Sources were centered in the slits by peaking up on 2MASS stars. 

The individual bed frames were processed by the S13.2.0 version of the SSC pipeline, 
which performs ramp fitting, dark current subtraction, droop and linearity corrections, flat- 
fielding, and wavelength and flux calibration^. Because the pipeline does not include back- 
ground subtraction, we have constructed, for the low-resolution long-slit spectra, a coadded 



2 See the IRS Data Handbook, 



http : //ssc . spitzer . caltech . edu/irs/dh 
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background frame from the bed observatio ns with the source in the opposing nod and off-order 



positions (see also IWeedman et aL||200£ ). Coadding was performed with the sigma-clipping 



option of the imcombine task in IRAEj. The inclusion of off-order and off-nod frames in 
the background image means that the total integration time on the background is three 
times that on the source, thus improving the signal-to-noise of the two-dimensional (2D) 
subtraction. 

We did not acquire a separate background spectrum for the high-resolution SH and LH 
spectra, and the slit size precludes using the same procedure as for the SL module. Therefore 
we subtracted the background from the SH and LH observations using the one- dimensional 
(ID) spectra as described below. The SH and LH 2D bed images were coadded in the same 
way as for the SL modules, and corrections for sporadic bad pixels and cosmic ray hits were 
carried out manually by inspection of the images at the separate nod positions. 

We extracted the source spectra with SPICE, the post-pipeline IRS package furnished 
by the SSC. Before extraction, all spectra at separated nod positions were cleaned with 



the IRSCLEAN algorithm (http://ssc.spitzer.caltech.edu/postbcd/irsclean.html). The au- 



tomatic point-source extraction window was used for all modules. Automatic extraction 
uses a variable-width extraction window which scales with wavelength in order to recover 
a constant fraction of the diffraction-limited instrument response. For the SL spectra, this 
gives a 4-pixel (7'.'2) length at 6 /im, and an 8-pixel length (14'/4) at 12 /im; the slit width 
is 3'.'6 for both SL modules. At high resolution, the SPICE extraction was performed over 
the entire slit (4'.'7xll / .'3 SH; ll'.'lx22'.'3 LH). Orders were spliced together by averaging, 



ignoring the noisy regions at the red end of each order (e.g., lArmus et al.l l2004j ) . Then the 
individual spectra were box-car smoothed to a resolution element, and clipped in order to 
eliminate any remaining spikes in the high-resolution data. Finally, the two spectra for each 
module (one for each nod position) were averaged. 

Background was subtracted from the ID SH spectra by defining a level from the overlap 
region with the SL spectra, from which we had subtracted a 2D background. Specifically, 
the difference between the 2D background-subtracted SL and SH spectra was minimized 
over the substantial region where they overlap in wavelength (~5/im). The spectral shape 
of the background was given by the model of Reach and coworkers^, and the multiplicative 
constant by the minimization process. Subsequently, the LH background was subtracted in 



3 IRAF is the Image Analysis and Reduction Facility made available to the astronomical community by 
the National Optical Astronomy Observatory, which is operated by AURA, Inc., under contract with the 
U.S. National Science Foundation. 



See http : / / ssc . spitzer . caltech. edu/documents/background 
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an analogous way, by minimizing the difference between the SH and LH spectra over their 
overlap region (<l/im). 

This "bootstrapping" procedure for background subtraction in the high-resolution mod- 
ules should in principle eliminate the need for scaling adjustments. However, its accuracy 
depends on the unresolved nature of the sources in the spectral apertures. We measured the 
source sizes by fitting Gaussians to the IRAC 7.9 /im and MIPS 24 /im images. The mean 
(and median) image FWHM is 3'.'7 at 7.9 /im, and 6'.'2 at 24 /im. The SL1 slit width is 3'.'6, 
so most of the sources are sufficiently compact to not significantly exceed the SL1 slit size. 
The resolution of MIPS24 is ~6", making the BCDs essentially unresolved at these longer 
wavelengths; in fact, the characteristic diffraction pattern is seen in all the MIPS24 images. 
The FWHM of the largest source at 24 /im (Tol 1924—416) is 10", so even in the most ex- 
treme case, the LH slit width of 11'! 1 completely encompasses the BCD light at 24 /im. Our 
assumption of point-like morphologies at Spitzer wavelengths thus appeared to be reasonably 
justified. Nevertheless, for a few (of the largest) sources, the level of the background needed 
to be adjusted slightly to provide a smooth transition among the modules. In any case, the 
24 /im flux as measured from MIPS is generally consistent with the continuum level derived 
for the IRS spectra (see below). 

We have obtained IRAC and MIPS photometry for all 23 BCDs in our sample. How- 
ever, we have IRS observations for only 22 BCDs. One BCD, SBS 0940+544, has an IRS 
observation in the Guaranteed Time program (PID85). Thus, we present here IRS results 
for only 22 objects, and will incorporate SBS 0940+544 into our sample in a future paper. 



2.2. IRAC and MIPS Imaging 



Although the focus of this paper is on the IRS spectra of our sample objects, we briefly 
describe their Spitzer photometric data because we need photometric quantities such as the 
total IR luminosity (Ltir) to properly interpret their spectral properties . A more d e tailed 



description of the photometry a nd reduction procedures can be found in iHunt et al.l ( 12006 



Haro3) and iThuan et al.l ( 120081 . Mrk996). The discussion of the photometric data of the 
whole sample is deferred to a future paper on SEDs. 

Briefly, we acquired IRAC images in two filters, 4.5 and 7.9 /im, each with two sets of 
four frames in the high dynamic range mode. The individual bed frames were processed 
with the S14.0.0 version of the SSC pipeline (see the IRAC Data Handbook^) . The bed 



5 Available from the SSC website [http : //ssc . spitzer . caltech . edu/irac/dh/ 
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frames (corrected for "banding" effects for bright sources) were coadde d using MOPEX, the 



imag e mosaicing and source-extraction package provided by the SSC ( IMakovoz &: Marleau 



20051 ). Pixels flagged by the masks were subsequently ignored. Additional inconsistent pixel 
values were removed by means of the MOPEX outlier rejection algorithms, in particular the 
dual-outlier technique, together with the multiframe algorithm. The frames were corrected 
for geometrical distortion and projected onto a "fiducial" (refined) coordinate system with 
pixel sizes of 1'.'20, roughly equivalent to the original pixels. Standard linear interpolation 
was used for the mosaics. The noise levels in our post-pipeline MOPEX IRAC mosaics are 
comparable to or lower than those in the SSC products. 

Our MIPS images were acquired in the Fixed Cluster- Offset mode in all three channels, 
with offsets of 12" in two additional pointings. The individual bed frames were processed by 
the the S14.4.0 version of the SSC pipeline, which converts the integration ramps inherent 
to the MIPS detecto rs into slopes, and corrects for temporal variations of the slope images 



(IGordon et al.ll2005l ). As for the IRAC images, we processed the dithered bed frames in 
the spatial domain with MOPEX. The DCE masks and the static masks were used to flag 
pixels which were subsequently ignored. The MOPEX outlier rejection was used to remove 
any additional spurious pixel values, with the dual-outlier and multiframe algorithm as for 
the IRAC frames. Geometrical distortion was corrected before projecting the frames onto a 
fiducial coordinate system with pixel sizes of l'/20 for MIPS24, roughly half of the original 
pixel size of 2'.'5. Pixel sizes of the final mosaics at 71 and 160 /zm are also approximately 
half the originals, i.e. 4'/95 at 71//m and 8'.'0 at 160 fim. Unlike the IRAC coadds, we 
incorporated the sigma-weighting algorithm because it gave less noisy MIPS mosaics than 
without. Standard linear interpolation was used in all cases. In all three channels, our 
MOPEX mosaics are superior to those provided by the automated post-pipeline reduction. 



2.3. IRAC and MIPS Photometry 



We have performed aperture photometry on the IRAC and MIPS images with the IRAF 
photometry package apphot, taking care to convert the MJy/sr surface brightness units of 
the images to integrated flux units. The background level was determined by averaging 
several adjacent empty sky regions. Fluxes were determined at radii where the growth curve 
has become asymptotically flat. In many cases, there are multiple sources within the IRS 
apertures, and we used subjectiv e judgment to match photometry with the spectroscopic 
slit. Following iDraine fc Lil (120071 ). these fluxes were used to derive the "color temperature" 
between 71 and 160 /im, (vL I/ (71))/(vL v (160)), and the total infrared flux (TIR): 



Ltir = 0.95 (uL v ) 7 j> + 1.15 (vL v ) u + (vL v ) n + (vL v ) 



Kin 
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Both quantities will be pertinent for the discussion of the spectra. 

In order to check our reduction and calibration procedures for both the IRS and MIPS, 
we have compared the MIPS 24 /zm total fluxes with the continuum levels in the IRS spectra. 
Such a comparison (see below) shows good general agreement, implying that both the IRS 
flux calibration (in the case of a point-source) and background subtraction and the derivation 
of the total MIPS flux have been done correctly. The only exception is UM311, which is in 
a crowded field (there are several point sources in close proximity to a spiral galaxy), and 
thus its photometry is probably contaminated by other sources. 

All targets were detected at 24 and 71 /zm. However, at 160 /mi, we obtained only upper 
limits (ULs) for 6 of them. As a consequence of a Spitzer queue problem, SBS 1030+583 has 
no MIPS24 observations. 



3. Spectral Analysis 



The flux emitted by aromatic (or PAH) features constitutes a large fraction of the 
IR energy budget in a typical star-forming galaxy spectrum. Ho wever, these features are 



known to decrease in intensity or even disappear at low metalli cities ( lEngelbracht et al.ll2005 



Madden et al.l l2006t IWu et al.l 120061 ; lEngelbracht et al.l 120081 ). Our sample contains many 



metal- deficient objects and is thus uniquely able to address the behavior of PAH emission 
at low metallicities. In order to directly compare the results for our low-metallicity sample 
with those for an existing well-defined set of more metal-rich galaxies, we wanted to analyze 
our IRS spe ctra in the same m anner as this latter set. This consideration dictated the use 
of PAHFIT (ISmith et al.l 120071 ) . an IDL pr ocedure which was dev eloped for and applied to 
the Spitzer Nearby Galaxy Survey (SINGS. iKennicutt et al.1120031 ). 



PAHFIT is particularly suited for separating emission lines and PAH features (e.g., 
the PAH blend at 12.6-12.7/im and the [Nell] line at 12.8 /mi), as well as measuring faint 
PAH features superimp osed on a strong continuum . Broad "plateaus" underlying the 8 /jjd. 
and 16-17 /im emission ( iPeeters et al.ll2002l . l2004bl ) and strong silicate absorption at 9.7/im 
can make the measurement of broad PAH bands exceedingly difficult in those wavelength 
regions. PAHFIT overcomes this problem by fitting simultaneously the spectral features, the 
underlying continuum, and the extinction, and gives quantitative estimates and uncertainties 
for all parameters. Drude profiles are fitted to the PAH features, and Gaussian profiles to 
the molecular hydrogen and fine-structure lines. In addition to PAHFITting the spectra, 
we also measured emission line fluxes over the entire IRS range by fitting Gaussian profiles. 
Our strategy for deriving quantitative measurements from the IRS spectra is described in 
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the following. 

3.1. PAHFIT 

Our spectra have a higher resolution than those normally fitted with PAHFIT, specif- 
ically for wavelengths >14.5 /im where we have no low-resolution spectra. Hence, after 
binning the spectra to a 0.03 fim resolution, and fitting the entire spectrum, we use only 
the spectral lines measured with PAHFIT having A < 14.5 /im. In contrast, we adopt all 
the dust features fitted by PAHFIT over the entire IRS range. We modified the PAHFIT 
IDL code to allow the central wavelengths and widths of the Drude profiles to vary, and 
experimented with different weighting schemes. We also tried several different combinations 
of dust continuum temperatures, but these did not influence the fit results. PAHFIT was 
run with three different options, including mixed extinction, screen extinction, and fixing 
extinction to zero. 

The nominally best fits (those with lowest xt) f° r H BCDs had non-zero extinction. 
However, in every case but Haro3 and Mrk996, the statistical significance of the fit was 
comparable to similar fits, but with zero extinction. Hence, when possible, we opted for the 
fit with no extinction [r(9.7) = 0]. 

The results are shown in Figure [TJ The top panels show the best-fit PAHFIT model as 
a red curve, and the bottom panels give the residuals. The reduced xl given in each panel 
in the figures is calculated from the residuals of the fit up to 14.5 /im, normalized by the 
mean spectral uncertainty over the same spectral range. The figures show that the residuals 
are generally quite well-behaved, and that xt * s generally sufficiently small to imply a good 
fit. Uncertainties on the integrated line fluxes were estimated in two ways. The standard 
deviation a per pixel of the continuum was measured individually in a relatively clean region 
of each SL, SH, and LH spectrum. A first estimate of the uncertainty in the line flux, 
assumed to be spread over 3 pixels, is taken to be 3 x 3a. This is compared with a second 
uncertainty estimate, the one given by PAHFIT for each spectral feature. The largest of 
the two estimates was then chosen to be the final uncertainty of the measurement. Table [2] 
reports all FS lines with A <14.5/im detected at the £ 3a level with PAHFIT, and Table [3] 
lists the dust features. 
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3.2. Gaussian Fitting of Emission Lines 

We fitted all possible molecular and FS emission lines with Gaussian profiles, taking 
the integrated fitted flux as the measure of the total flux in the line. When necessary, we 
used multiple Gaussian profiles to deblend adjacent lines (e.g., [Oiv] and [Fell] at ~26//m). 
Linear continua were fitted simultaneously, taking care to define the continuum regions well 
beyond the Gaussian wings. 

In order for a spectral feature to be identified with a known emission line, its rest- 
frame wavelength shift from the laboratory wavelength had to to be <0.065 /im (for each 
object, we correct the observed spectrum to its rest frame by using the optical redshift 
listed in Table [T|) . Each emission line was measured several times, varying the continuum 
positions and blending parameters, to obtain an additional estimate of the uncertainty. The 
uncertainty on the integrated flux was then taken to be the larger of the two following 
estimates: the standard deviation of the ensemble of different measurements, or 9 times the 
standard deviation a of the continuum fluctuations, as described above (assuming the line is 
spread over 3 pixels). To guard against false detections from bad pixels or spectral glitches, 
we also required that the average of the two nod positions have a signal-to- noise >3. Table 
H] reports the long- wavelength FS fluxes we obtain from these measurements. 

For the objects where [O iv] or [Fe il] was detected, blown-up parts of the spectra showing 
the 26/im region are given in Figure [2j Although we detect [O iv] at 3a (2a) in 7 (9) galaxies, 
[Fell] is detected only in 2 (5 at 2a). This result will be discussed in more detail in §4.41 
The extremely high-excitation lines [Nev] (ionization potential of 97.1 eV) at 14 and 24 /im 
are not detected in our data for any BCD. The MIPS 24 /im continuum level is also shown 
in Fig. [2] by an open circle. As stated in §2.3} the IRS spectra and MIPS 24 /im total fluxes 
are consistent, except for UM311 which is in a crowded field. 

In addition to the definite presence of H 2 , a few spectra shown in Fig. [3] also tentatively 
suggest the presence of two other molecules, OH (~28.9, 30.3, 30.7 /xm) and H 2 (~29.8 
and 29.9 /im). Table gives the H 2 line fluxes (including those measured by PAHFIT for 
A < 14.5 /xm, see above), and Table [6] reports the other molecules. These detections will be 
discussed further in §6j 

3.3. Comparison of Fitting Techniques 

To judge the reliability of our line flux measurements, we have compared the fluxes of 
those emission lines that fall in the short- wavelength region of the spectra (A < 14.3 /im), and 
that have been measured in two different ways: first by PAHFIT and second by individual 
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Gaussian fits. For strong lines with fluxes ~ 2 x 10~ 17 Wm -2 , the difference between these 
two methods is 4%; for faint lines with fluxes ~ 5 x 10~ 18 Wm~ 2 , the difference can be 
as large as 25%. The larger discrepancy for the fainter lines can be attributed to the choice 
of different continuum levels in the two methods. This is the main reason why we preferred 
to use PAHFIT for the short-wavelength region. However, the overall agreement between 
the two methods is sufficiently good to conclude that we can reliably compare short- and 
long-wavelength line fluxes, and that we do not introduce systematic differences in measuring 
them with slightly different methods. 

Although PAHFIT corrects for dust extinction, we have not applied an analagous cor- 
rection to the line fluxes for A >14.5 /im. Such corrections would be applicable only to two 
sources, Haro3 and Mrk996. Haro3, with r(9.7) = 0.39, would have a correction 25% for 
wavelengths longer than 15 /im; Mrk996 would have a 2% correction. Hence, because the 
correction for both objects are comparable to or smaller than uncertainties in the line fluxes, 
we did not apply them to the long-wavelength measurements. 



4. The Ionized Component of the ISM 

The spectral region covered by Spitzer/IRS provides a wealth of diagnostics to probe the 
physics of the ISM at low metallicity. The FS lines allow us to assess the physical conditions 
of the ISM, including the spectral shape or hardness of the ISRF and the density of the ionized 
gas. These quantities affect the dust properties by causing changes in the temperature of 
the "classical" grains, and by altering the stochastic emission from the aromatic features. 
They constrain the molecular content through dissociation and excitation processes. 



4.1. The Hardness of the Interstellar Radiation Field 



The IR FS lines span a wide range of ionization potentials and are thus sensitive diag- 
nostics of the ionized ISM in BCD Hll regions. We focus first on diagnostics of the ISRF 
hardness. In their discussion, we will attempt to distinguish among the interdependent 
effects of ISRF hardness, ISRF intensity, and nebular metallicity. 

Traditionally, flux ratios of high-ionization species have been used to measure the ISRF 
hardness, or equivalently the UV slope of its spectrum, which is sensitive to the effective 
temperatures of the ionizing sta rs. In the IR range, the [Ne ill] /[Nell] an d [Siv]/[Siii] flux 



ratios are usually adopted (e.g., iThornley et al.ll2000t iGiveon et al.ll2002l ). Such flux ratios 
also depend on the ionization parameter U, defined as U = Q / AnR*n e c, where Q is the 
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number of ionizing photons, n e is the ele ctron density, and i?* is the radius of the star- 



forming region (e.g., iThornley et al.ll2000l ). We first examine whether these ratios depend 
on metallicity as shown in Fig. |H The top panel shows that within our sample there may be 
a very slight dependence of the [Siv]/[Sm] ratio on metallicity, although with a very large 
scatter (there are no SINGS data for the variation of this flux ratio with abundance). The 
bottom panel shows that there is very little correlation of the [Ne ill] /[Nell] flux ratio with 
oxygen abundance within our sample alone. However, when the SINGS data are included, it 
is clear that the low-met all icity BCDs show considerably higher neon ratios than the more 



metal-rich SINGS galaxies (IDale et al.ll2009l ). The SINGS galaxies span an abundance range 
of roughly a factor of 10, from 12+log(0/H) ~8.4 to 9.4, and our sample fills in the decade 
below, down to 12+log(0/H) ~7.5. Considering the two samples together, the bottom panel 
shows an increase of the [Ne in] /[Neil] flux ratio by a factor of ~10 from 12+log(0/H) ~9.4 
down to ~8.2. Then the curve flattens out into a "plateau" for lower abundances, with 
[Ne m]/[Ne il] ~5. Again, the scatter of the [Ne m]/[Ne il] flux ratio at a given oxygen 
abundance is large, indicating that metallicity cannot be the main parameter controlling 



(Wu et al. 


2006; 


Hao et al. 


2009) 



Our sample was defined on the basis of the optical Hell/H/3 ratio which gives an inde- 
pendent estimate of the ISRF hardness. However, within our sample, there is no correlation 
between the [Ne ill] /[Nell] and [Siv]/[Siii] ratios and the Hell/H/3 ratio as shown in Fig. El 
A possible reason for this lack of correlation is that the IR and optical ratios do not probe 
the same excitation energy. The ionization potentials for [Ne ill] and [Ne il] are respectively 
41.0 eV and 21.6 eV, while they are respectively 34.8 eV and 23.3 eV for [Siv] and [S ill] . 
These ionization potentials are well below that of He + which is 54.4 eV. Thus, the ISRF UV 
energies probed by the neon and sulfur lines are too soft compared to those probed by the 
Hell ratio. 

We have therefore examined the relation between the [Olv]/[Sill] and Hell/H/3 ratios; 
the ionization potential of [O iv] is 54.9 eV, similar to that of He + (54.4 eV), while that of 
[Sill] is only 8.2 eV. Since both silicon and oxygen are products of core-collapse supernovae 
(SNe), the [Olv]/[Sill] ratio should not be subject to abundance anomalies. The left panel 
of Fig. [HI shows a weak correlation between the two ratios for the BCDs in our sample: a 
high [Olv]/[Sill] ratio is associated with a higher Hell/H/3, at a confidence level of ~93% 
(one-tailed). Higher [Olv]/[Sill] is associated with lower nebular O/H abundance (right 
panel of Fig. [5]), at a ~98% confidenc e level (one-tailed). T his last correlation is similar 
to the one between Hell/H/3 and O/H (IThuan fc Izotovl 120051 ) . Evidently, lower metallicity 
galaxies tend to have generally harder ionizing radiation. 
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Although the BCDs by themselves show little correlation between the [Oiv]/[Sin] ratio 
and the [Ne in] /[Nell] and [Siv]/[Siii] ratios, when SINGS galaxies are included, there is 
a definite correlation. Fig. [7] illustrates this, and suggests that in star-forming galaxies the 
presence of hard radiation as probed by the [Olv]/[Sill] ratio implies also the presence of 
less hard radiation as probed by the [Ne ill] /[Nell] and [Slv]/[Sm] ratios. 



4.2. The Intensity of the Interstellar Radiation Field 



We now explore diagnostics for the ISRF intensity, as distinct from the ISRF hardness, 
since it does not depend on the spectral shape of the ionizing radiation. The ISRF intensity 
is related to the ionization parameter which measures the number of ionizing photons per 
hydrogen atom. MIPS observations are relevant here, because the ratio of 71 to 160 /im 
fluxes is sensitive to the temperatu re of the large ("cl assical") grains, and can thus act as 
an indicator of the ISRF intensity (jDraine fc Lil 120071 ). Essentially, we will use the MIPS 
luminosity ratio ("color temperature") (z/Lj,(71))/(z/Lj,(160)) as a proxy for the ionization 
parameter. The upper panel of Fig. [8] shows the relation between this "color temperature" 
(i/L„(71))/(i/Lj,(160)) and the quantity {vL v )24,/{{vL v )yx-\- (vL„)i6o); the latter represents 
the relative contribution of the 24 ju m emission to the lo ng-wavelength component of I/tir, 
the total IR luminosity as defined by lDraine &: Lil (120071 . see §2.3p . The lower panel of Fig. 



M shows the variation of the same 71 to 160 /im luminosity ratio with oxyge n abun dance. 
SINGS galaxies are also shown in Fig. [8] using data taken from iDale et al.l (120071 ). The 
correlations in both panels are highly significant even within the BCD sample alone: the 
correlation in the top panel is at a £ 99% confidence level (one-tailed), while the bottom 
one is at a £ 99.9% confidence level (one-tailed). This indicates that the 24/zm fraction 
of (i/L w (71)) + (z/Lv(160)) increases with increasing large-grain temperature (and decreasing 
oxygen abundance). Hence, the relative strength of the 24 /im luminosity is a good tracer of 
the intensity of the ISRF. As dust gets warmer, more large -grain emission is shifted toward 



shorter wavelengths, resulting in increased 24/zm flux (e.g.. lDale et al.ll2005l ) 



Apparently, at low metallicities, 12+log(0/H) 7.8, large grains can be heated to ex- 
treme (z/Lv(71))/(zaLv(160)) ratios, reaching values of 10-30. The effect could be even more 
extreme since the lowest metallicity BCDs in our sample are not plotted; the MIPS 160 fim 
emission for th e four galaxies with 1 2+log(Q/H) £ 7.6 is not detected. Compared to the 
SINGS sample JPale et alJ2005ll2007h . the ISRF intensity as measured by (i/L I/ (71))/(z/L„(160)) 
is 5 times greater on average in the BCDs than in even the most intense SINGS galaxjQ. 



3 To do the comparison, we have converted the flux ratios used bv IDale et al.1 (|2005l ) to the {vL v ) ratios 
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Using the 71/160/xm luminosity ratio ((vL„(71))/(i/L,,(160))) as a tracer of the ISRF 
intensity, we next investigate whether any of the FS line ratios, indicators of the ISRF 
hardness as shown above, are correlated with it. Figure M shows the [Ne ill] /[Nell] and 
[Oiv]/[Sin] line rati os as a function of the 71/160£tm r a tio; S INGS galaxies are also plotted 



combining data from Dale et all ( 120071 ) and iDale et all (120091 ). The [Ne ill] /[Nell] flux ratio 
shows no correlation, although there is a weak one in the BCD sample alone (at the £ 97% 
confidence level) for [O iv]/[Sin]; the [Olv]/[Sill] ratio is larger for warmer large grains. 
In the SINGS sample, most of the highest [Ne in] /[Nell] and [Olv]/[Sill] flux ratios are 
associated with AGN, so any correlation could be masked by AGN high-excitation emission- 
line ratios in galaxies with a relatively low-intensity ISRF. In any case, we conclude that 
the hardness of the ISRF, as measured by [Olv]/[Sill] is weakly related to its intensity, as 
traced by the 71/160/zm luminosity ratio. In lower metallicity BCDs, the ISRF tends to be 
both harder and more intense. 

To summarize, in our low-metallicity BCD sample, the [Olv]/[Sill] ratio is correlated 
with Hell/H/3, 12+log(0/H), and the large-grain dust temperature (71/160/xm ratio). This 
implies that star-forming galaxies with a lower nebular metallicity as compared to more 
metal-rich galaxies, generally possess an ISRF which is both harder, as indicated by higher 
[Olv]/[Sill] or Hell/H/3 flux ratios, and more intense as indicated by a larger 71/160/xm 
luminosity ratio. The relatively large scatter in some of these relations, and the "plateau- 
like" behavior of others, implies that these quantities do not depend on only one parameter, 
but that many parameters are involved. The ISRFs in the SINGS galaxies are of considerably 
lower intensity than those generally found in BCDs. This is probably because the star- 
forming regions are more extended, less extreme, and have a lower ionization parameter Q. 
Such factors are probably as important as ISRF hardness and metallicity. 



4.3. Electron densities 



We have derived electr on densities n e in the ionized gas from the ratio of the temperature- 
insensitive [S in] IR lines ( iDraind 120091 ). Fig. [TU1 shows these densities as compared with 
those derived using the optical [Sll] A6717/A6731 ratio. If the IR emission were to originate 
in regions suffering from higher extinction than the optical regions, we might expect the IR 
n e to be higher than those derived from the optical lines. Fig. [TD] shows that this is true 
for some BCDs, but not all. On the contrary, a few galaxies show higher optical than IR 
densities. The reason for this behavior is unclear; beam dilution could play a role since the 



used here. 
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IRS aperture is several times larger than the optical one and thus includes lower density 
outer regions. 

In the optical, high-ionization radiation responsible for ions such as [Nev] is found to 
be associated with extr emely massive compact super star clusters and high optical densities 
( jThuan fe Izotovll2005l ). The idea that compact star- forming regions we re associated with 



higher electron densities n e was first proposed on theoretical grounds by iHirashita fe Hunt 



( 120041 ) . The size-density correlation in Hn regions in which smaller region s are denser is ob- 
serve d over six orders of magnitude from Galactic objects to distant BCDs (jHunt fe Hirashita 
2009 J). If fast radiativ e shocks were responsible for the ionizing radiation as proposed by 
Thuan fe Izotovl ( 120051 ). then we would expect such processes to be more efficient at high 
densities. However, in the current BCD sample, only two (three) of 7 (9) galaxies with a 3cr 
(2<r) [O iv] detection has a n e above the low-density limit of ~30 cm~ 3 . The electron densities 
as derived from the [S ill] IR lines are apparently not correlated with the hardness of the 
ISRF, as measured by [O iv]/[Sin]. One reason for this non-correlation may be that [Oiv] 
with an ionization potential of 54.9 eV is not probing sufficiently high ionization energies 
(the [Ne v] line has an ionization potential of 97.1 eV). Another reason may be that the [S ill] 
IR lines arise from a region that is less dense than the [O iv] region. More work on a larger 
sample is needed to explore the connection between density, geometry, and the ISRF. 



4.4. [O IV] and [Fe il] in BCDs 



The two FS lines [O iv] and [Fell] lie at adjacent wavelengths (AA =0.098/im) near 
26.0 /im. They are, however, sufficiently separated to be resolved by LH IRS spectra. While 
[O iv] traces high-excitation gas, [Fe il] has an ionization potential of 7.9 e V, more similar to 



[Sill] (see ^4.51). Both l ines may be excit ed by photoionization and shocks ( ILutz et al.lll998 



O'Halloran et al.ll2006l ; lAUen et al.ll2008l ). By fitting the two lines with Gaussian profiles, we 
are able to deblend them in our spectra. We detect [Oiv] at 3a in 7/22 observed galaxies, 
but [Fell] at 3a in only Mrk5 and Mrk996 (see Fig. [2]). At low metallicity, [Oiv] emission 
is almost 4 times as common as [Fell] (32% vs. 9%). 



This conclusion contrasts with the findings of lO'Halloran et al.l (120061 120081 ). Using 
Spitzer archival IRS data, they have published 26.0 /im [Fe il] fluxes for a sample of starbursts, 
including the 22 BCDs reported here. While the uncertainties in the spectra estim a ted b y 
us (see §|3]) are generally more conservative than those reported by lO'Halloran et al.l ( 120081 ). 
these authors find [Fe il] emission in all objects, but report no [O iv] emission. The [Fe il] 
emission fluxes they find are also larger than ours, suggesting contamination from the [O iv] 
emission. Figure [2] shows the identifications of the two lines; it is evident that the [Fe il] flux 



- 16 - 



can be securely measure d only when [Oiv] is taken into account. No spectra are shown by 



O'Halloran et al.l (120081 ) . so we cannot directly compare their line identifications with ours 



to resolve the source of discrepancy. 



4.5. [Sill] Emission 

The [Si il] line is detected in more than half the BCDs in our sample. At an ionization 
potential of only 8.2 eV, [Sill] traces the neutral gas, rather than the more highly excited 
ionized gas traced by the neon, sulfur, and oxygen lines. Figure [UJ shows the [Si n] emission 
normalized to the TIR luminosity. The mean [Sill]/TIR ratio is more than a factor of two 
lower than that of the SINGS galaxies, 8 x 10 -4 versus 2 x 1CT 3 (compare the dotted and 
dashed horizontal lines in Fig. [TT]) . This slight difference is probably due to the lower metal 
abundance of the BCDs, and the consequently less efficient cooling of the gas through FS 
lines. It could also be due to a more intense ISRF, and thus a lower Si + abundance. In any 
case, the ratio of [Sill] emission to Ltir is relatively constant with luminosity (excluding 
UM311 as discussed previously). These conside rations suggest that the [Sill] line is an 



efficient coolant, comparable to the [Cn] line (e.g.. iMalhotra et al.lll997l ). and could be used 



as a proxy for Ltir to within a factor of two, even at low metallicities. 



Aromatic Features 



The wealth of spectral data obtained by ISO combined with much theoretical work 
and laboratory measurements has shown convincingly that aromatic features are associ- 



ated with stochastic or single- photon emission from large molecu 



matic hydrocarbon s or PAHs (|Joblin et al 



2008, 



Peeters et al. 2002; van Diedenhoven et al. 



1995 



2004 : 



Verstraete et al. 



Peeters et al- 



es cal l ed polycyclic aro- 



2001 



2004a : 



Honv et al 



2001 



Bauschlicher et al. 



2009h. PAH emission dominat es the mid-infrared spectra of star-forming galaxies 



( iBrandl et al.ll2006l ; ISmith et al.l 120071 ). and contributes importantly to their IR energy bud- 
get. Through the photoelectric effect, PAHs heat the gas in the ISM (e.g.. lHollenbach fc Tielens 
19971 ). thereby reinforcing the coupling of the gas and dust components. Because they are so 
dominant in the MIR spectral range, PAHs have been used to identify galaxies with intense 
star formation at redshifts of ~ 2 — 3, and to estimate their bolometric luminosities (e.g., 



Houck et al.l 120071 : IWeedman k Houckll2008l : iDev et al.l 120081 ) 



In most star-forming galaxies, PAH emission shows a surprisingly narrow range of 
properties, well characterized by a "standard set" of features with fixed wavelengths and 
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FWHMs (ISmith et all 120071 ). However, first ISO and lat er Spitzer have shown that metal- 



Eneelbracht et al. 


2005: 


Madden et al. 


2006 



though it has generally been concluded that the reason for this deficiency is low metallicity, 
it is not yet clear whether it is the metallicity directly, or rather the indirect effects of a low 
metallicity environment. In the first case, there would be a consequent lack of raw material 
(carbon and nitrogen) with which to form PAHs; in the second, the increased ISRF hardness 
or intensity would lead to their destruction. Here, we examine in detail the PAH emission 
in our sample to examine this question, and show for the first time that some PAHs can 
survive even in extremely metal-poor environments. 



5.1. PAH Properties at Low Metallicity 

The 7.7 jum blend is the most common PAH feature in our sample, being detected in 
15 of 22 objects. The 7.7 /xm feature is also the strongest one, comprising on average 49% 
of the total PAH power. The remaining PAHs are significantly weaker, and less frequently 
detected: 13 BCDs show the 11.3 /xm blend, 9 the 8.6 /xm feature, and 7 have a detection 
of the 6.2 /xm PAH. No 17 /xm PAH features were detected, and only Mrkl329 showed a 
significant detection at 16.4 /xm. The 6.2, 7.7, 8.6, 11.3, and 12.6 /xm features by themselves 
constitute ~72% of the total PAH power in BCDs (see Table El); the SIN GS sample has 



about ~85% of the total PAH power in these features (ISmith et all 120071 ). Similarly to 
the SINGS galaxies, other weaker features (not listed in Table |3]) which contribute to the 
PAH luminosity include the 5.7, 6.6, 6.7, 8.3, 10.6, and 10.8 /xm features. In our sample, 
these features comprise roughly 28% of the PAH luminosity, but in SINGS galaxies, only 
15%. Part of the reason for this difference may be the difficulty in detecting long-wavelength 
PAHs in BCDs; only 1 BCD has the 16.4 /xm feature detected, while in SINGS galaxies the 
16-17/xm blend makes up ~6% of the total PAH power. 

The fractional power of the four strongest features relative to the total PAH luminosity 
is illustrated in Figure [12j The horizontal lines in each panel correspond to the SINGS 
medians (dashed) and the BCD means (dotted). The BCD means are calculated taking into 
account all galaxies with 7.7 /xm detections; thus they are a sort of weighted average which 
considers frequency of detection together with intensity. This is the reason that the mean 
for the 6.2 /xm PAH lies below most of the BCD data points: that feature was detected 
only in 7 galaxies, while the other features in Fig. [12] were detected with a frequency more 
similar to the 7.7/xm PAH. The relative PAH fractions for the BCD and SINGS 6.2 and 
7.7 /xm features are virtually indistinguishable: 0.10 vs. 0.11 for 6.2/xm, and 0.49 vs. 0.42 
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for 7.7 /im (ISmith et al.ll2007l ). However, for the longer wavelength PAHs at 8.6 and 11.3 /im, 
the BCD fractions are > 40% larger: 0.10 vs. 0.07 (8.6) and 0.17 vs. 0.12 (11.3). Although 
there is considerable scatter, the BCD mean PAH relative fractions for these features tend 
to exceed even the 10 to 90th percentile spreads of the SINGS galaxies. 



Figure [13] shows the mean Drude profile widths and central wavelengths of six aromatic 
features detected at 3cr in our sample; the numbers of BCDs with each feature are given 
in parentheses. The analogous quantities for the SINGS galaxies are also plotted. In no 
case, were the PAHFIT profiles fixed to the SINGS galaxy parameters able to fit the BCD 
data; every BCD was better fitted by allowing the Drude profile widths (FWHMs) and 
central wavelengths to vary. In some cases, the differences between the xl with the SINGS 
"standard" profile parameters and the best-fit parameters were small, £ 15%; but in others, 
the improvement of \\ was almost a factor of 2. The only systematic difference between 
the two samples is the narrower profile width for the 8.6 /im feature; only the very broadest 
BCD profiles are as broad as those in more metal-rich systems. 

The larger relative intensities of the 8.6 and 11.3 /im feature in BCDs and the nar- 
row profile width of the 8.6 /im PAH could be due to a differen t size distribution of PAH 
populations at low metal abundances. iBauschlicher et all (120081 ) find that the 8.6 /im band 
arises from large PAHs, with A^ m i n £ 100 carbon atoms. Hence, they suggest that the rel- 
ative intensity of the 8.6 /im band can be taken as an indicator of the relative amounts of 
large and small PAHs in a given population. The large 8.6 /im intensity relative to total 
PAH power in BCDs could be a signature of fewer small PAHs (or more larger ones) at 
low metallicity. The relative lack of small PAHs could also explain the low detection rate 
of the 6.2 /im feature, sinc e most of its intensit y comes from PAHs with less than 100 C 
atoms ( Schutte et al. 1993 ; Hudgins et aL 12005). The siz e of a PAH molecule grows with 



the number of carbon atoms in it (e.g., iDraine fc Lil 120071 ) . For a given maximum number 

nber N m i n theoretically gives narrower profiles 



(Joblin et al. 


1995; 


mx , „ — .,„, ..„„.„ 
Verstraete et al. 


2001) 



profiles by symmetric Drude profiles in PAHFIT is a simplification, the narrower width of 
the 8.6 /im feature could be a signature of relatively larger PAHs at low metallicity. The 

engths than 8.6 /im 



difference disappears at 7.6 /im and is less pronounced at longer wave 
because both small and large PAH sizes contribute to these bandjll (Schutte et al. 1993 



Peeters et al. 2002 ; van Diedenhoven et al. 2004 ; Bauschlicher et al. 20091 ). 



Another indication that low-metallicity BCDs may be lacking the smallest size PAHs 



7 Thc 7.8 /im band may be dominated by larger PAHs only (j Bauschlicher et al, 
present in only two of the BCDs in our sample, Haro3 and IIZw70. 



2009); but this band is 
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comes from a correlation analysis. In a detailed study of Galactic Hll re gions, young stellar 
objects, reflection and planetary nebulae (RNe, PNe), and evolved stars. iHony et al.l (120011 ) 
found only a few correlations among PAH features. One of these is between 6.2 and the 
12.7 fim emission, which arises from the CC stretch mode in ionized PAHs, in co ntrast to the 
11.3 a m band which is at tributed to the neutral CH out-of-plane bending mode (IHony et al. 



2001 



Peeters et al.ll2002l ). The comparison between the 6.2/11.3 and 12.7/11.3 flux ratios 
is shown for our sample of BCDs in Fig. [T3J we show galaxies with £ 3a detections in at 



least two features. With the exception of Mrkl315, galaxies with 6.2 /im PA H detections 



and_12+log(0/H)>8.1 (filled circles) are similar to the Hll regions studied by IHony et al. 



( 1200 ll ). On the other hand, BCDs with lower O/H occupy a region of the plot with small 
6.2/11.3 ratios. Although the statistics are poor, small 6.2/11.3 are found for most of the 
lowest-met a l licit v BCDs. They lie in the same region as the RNe and PNe studied by 



Hony et al.l (120011 ) . Because the 11.3/im feature arises from non-adjacent or "solo" CH 
groups, dominant 11.3/im emi ssion implies large PAH complexes, with £ 100-200 carbon 



atoms and long straight edges (IHony et al.l 12001c iBauschlicher et al.ll2008l 120091 ) . Moreover 
fewer small PAHs reduce the intensity of the 6.2 /im e mission, as the bulk of this feature 



come s from PAHs containing less than 100 C atoms (jSchutte et al.l 11993c iHudgins et al. 



20051 ). Therefore, small 6.2/11.3/im band ratios at low metallicites could imply an overall 
deficit of small PAHs. 

A third indication that there is a deficit of small PAHs in a metal-poor ISM is provided 
by the intensity of the 8 .6 /im feature relative to the 7.7/im blend. As mentioned above, 
Bauschlicher et al.l (120081 ) find that 8.6 /im band arises from large PAHs, with iV min £ 100 



carbon atoms. O n the other hand, the 7.7/im feature contains emission from both sma ll and 
large PAH sizes (jSchutte et al.l Il993l ; iPeeters et al.l |2002| ; Ivan Diedenhoven et al.l 12004 ) . For 
the 9 BCDs with the 7.7 and 8.6 /im PAH features detected at £ 3a, the mean 8.6/7.7 flux 
ratio is 0.48, and the median is 0.33. This ratio is at least double that in the mo re metal-rich 
SING S galaxies, with 8.6/7.7 ~0.18 in the mean, and a range of 0.11—0.21 (jSmith et al. 



20071 ); 7 of 9 BCDs are outside the upper 90th percentile limit of 0.21. Again, although 
the number statistics are small, the large 8.6/7.7 /im flux ratios of the PAHs in these low- 
metallicity BCDs appear to be dominated by the largest complexes. 

In summary, there are three different lines of evidence which tentatively suggest that 
the PAH populations in low-metallicity BCDs are deficient in the smallest sizes (iV m i n £ 50 
C atoms). The IRS spectra show flux ratios typical of the largest PAHs modeled so far, with 
Amm ~ 100 C atoms. Apparently, the ISM environment in the BCD Hll regions is not pro- 
pitious to the existence of the smallest PAHs and has led to their destruction. Perhaps only 
the PAH population with the largest complexes can survive the extreme physical conditions 
in the low-metallicity Hll regions in the BCDs. 
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5.2. PAHs and the Energy Budget 

PAHs emit about 10% of the total infrared power in typical s olar metallicity star -forming: 



galaxies (12+log(0/H) ~ 8.7-9.0), with a maximum of ~20% flSmith et al.ll2007f ). In our 
low-metallicity BCD sample, with oxygen abundances 12+log(0/H) ranging from 7.4 to 8.3, 
the fraction of PAH emission ranges from ~0.1% to ~1.6%, with a median (and mean) of 
~O.5%0. Figure [T5l shows how the sum of all ^ 3cr PAH emission features normalized to 
the total IR luminosity Ltir, S(PAH)/TIR, varies as a function of Ltir for our BCD (filled 



circle s) and SINGS galaxies (non-AGN open squares, AGN open triangles) (ISmith et al. 



20071 ). For the SINGS galaxies, Ltir has been calculated within the spectral extraction 
aperture, while for the BCDs, Ltir is taken to be the total luminosity as described in §2.31 
Because most of the BCDs are point-like at 24 jum and are generally unresolved with MIPS, 
this is usually comparable to Ltir wit hin the spectral ap erture (see Fig. [2]). Were the source 



sizes larger at 160 /mi (e.g., Haro3: iHunt et al.l 120061 ). Ltir could be overestimated, and 



the PAH fractions somewhat higher. The infrared luminosities Ltir of the SINGS sample^] 
vary by a range of more than three orders of magnitude, while the BCD luminosities occupy 
the middle portion of this range. Figure [15] shows clearly the segregation of the SINGS 
points, with high E(PAH)/TIR, and the BCD points with low S(PAH)/TIR. The AGN, 
with intermediate S(PAH)/TIR, lie between the metal-rich star-forming SINGS galaxies 
and BCDs. Figure [151 also shows that there is no trend of the PAH fraction S(PAH)/TIR 
with Ltir m either sample. 

We might expect the fraction of IR luminosity in PAHs to vary with metallicity. We 
explore this possibility in the left panel of Fig. [161 where E(PAH)/TIR is plotted against 
oxygen abundance. Although there is no correlation of PAH fraction with metallicity within 
the BCD sample, the inclusion of SINGS galaxies does show a trend of increasing PAH 
fraction with increasing metallicity. We next investigate other physical factors which may 
control the PAH fraction, such as the ISRF hardness. The right panel of Fig. [16] shows 
E(PAH)/TIR as a function of the IR neon line ratio, [Ne in] /[Neil] (only those BCDs with 
at least a 3a detection in the neon lines are shown). As discussed in §4] the [Ne ill] /[Nell] 
ratio is a good tracer of ISRF hardness, up to an intermediate energy range. There is 
a weak correlation between S(PAH)/TIR and [Ne ill] /[Nell], significant at the 95% (one- 



8 Because of crowding and source structure, the TIR luminosity in UM311 may be overestimated; never- 
theless even without UM311, the mean is only slightly larger, 0.57% vs. 0.54%. 

9 The normalizing factor, Ltir, is calculated by different groups with slightly different algorithms; however, 
the general consistency of the data suggests that our conclusions do not depend on the exact method used 
to derive Ltir- 
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tailed) confidence level. Thus, over the range in UV energies £ 41 eV, the fraction of Ltir 
emerging as PAH emission may depend weakly on the hardness of the ISRF. This suggests 
that the PAH emission fraction in BCDs is not directly controlled by metallicity, but by 
the hardness of the radiation field. The latter could also be responsible for the destruction 
of the smallest PAH particles, as discussed in the previous section. We have checked for a 
correlation of S(PAH)/TIR with ISRF intensity as measured by (uL u (71))/{uL ll (160)), but 
found none. 

To conclude our discussion of PAHs and the IR energy budget, we examine the di- 
agnostic proposed by iHouck et al. J2007h : the quantity vL v (7.7[im), where L u (7.7 fim) is 
determined from the flux density at the peak of the 7.7+7.8 /im blend, with no correction 
for the underlying continuum. Because it arises from a range of sizes, the 7.7/xm band is 
arguably the most reliable PAH diagnostic, as it avoids preferential selection of either small 
or large PAH populations (ISchutte et all Il993l ; iPeeters et al.l l2002h . Figure [T7l compares 



^L,,(7.7/im)/Z/TTR for o ur BCDs with that for a sample of local starbursts (IBrandl et al. 



20061 ; IHouck et al.l 120071 ). Considering only those BCDs with a 160 /im detection, the BCDs 
have a mean vL v {7 .7 [im) / L^yr roughly four times lower than that of the starbursts. The 
means of the two samples are shown by horizontal dotted lines. Part of this difference may 
be due to a luminosity effect since the BCDs are more than 100 times less luminous in the 
mean than the starburst galaxies (see for example the SINGS galaxies in Fig. [T5|) . It may 
also result from variations of physical conditions in the ISM, since the BCDs with the high- 
est uL u (7.7/im)/L T i R clearly overlap with the range observed in local starbursts that are 
however much more luminous in Ltir. 



5.3. PAH Diagnostics for AGN/Starbursts Revisited 

Several diagnostics based on IR spectra have been proposed to separate AGN from star- 
forming regions in galactic nuclei. Previous work has also focused on the comparison of the 
flux ratios of PAH features and emission line ratios from FS transitions. In this section, we 
examine our BCD sample in the context of such diagnostics and assess their usefulness in 
the low metallicity regime. 



5.3.1. 7.7/11.3 nm Band Ratios 



We saw in the previous section that the total PAH fraction E( PAH)/Lttr i s corre lated 
with the IR neon ratio [Ne ill] /[Nell] (Fig. [TB]). Here, following ISmith et al.l ( 120071 ) . we 
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examine the variations of the 7.7/11.3 /im band ratios with [Ne ill] /[Nell] for both the BCD 
and SINGS samples as shown in Fig. [THJ Among the SINGS galaxies, there is a weak 
correlation, mostly for AGN, in the sense that a s maller 7.7/11.3 ratio is associated with 



a larger [Nelll]/[Nell]ratio (e.g., ISmith et al.l 120071 ). The BCDs extend the correlation to 



considerably higher neon ratios (by one order of magnitude), that is to harder ISRFs. The 
plot shows that the radiation fields in the low metallicity BCDs can be as extreme as those 
in the SINGS AGN. Such hard, intense ISRFs could suppress the ionized 7.7 /zm feature, 
relative to the neutral 11.3 /im feature. This may be due to the dominant contribution of 
large PAHs to the 11.3 /im band, as compared to the mixed size distribution thought to be 
responsible for the 7.7 //m band; the small PAH molecules contributing to the 7.7 /im emission 
could be destroyed while the large PAH molecules responsible for the 11.3 /im emission are 
not (see §5.ip . 



5.3.2. AGN / Starburst Diagnostics 



Because a hard, intense radiation field changes the PAH properties, a deficit of PAH 
emission, combined with FS line indicators of excitation, has often been used as a AGN 
diagnostic in solar metallicity star-forming galaxies. However, we have seen above that the 
radiation fields in low metallicity BCDs can be just as intense as those in AGN (see Fig. [IS]) , 
possibly destroying the smallest PAHs. Thus a PAH emission deficit may not necessarily be 
a reliable diagnostic for AGN. 



We compare our BCD sample with AGN ta ken from IGenzel et al.l (119981 ). using the 



;aurent et al. 



Genzel et al. 



(120001) . Figure QJ shows 
1998) as a function of the 



diagnostics proposed by the latter authors and 
the [Olv]/[Nell] ratio (or 1.7x33.5/im [Sill], see 
"strength" of the 7.7/ahi PAH band. [Oiv]/[Neii] is another measure of ISRF hardness (see 
§|4] where we used [O lv]/[Si n]). The "strength" of the 7.7/im PAH feature was determined 
by estimating the continuum level between 5 an d 14 /an and d ividin g the PAH flux by this 
level; this method is the same as that used by IGenzel et al.l ( 119981 ). so we preferred it to 
the simpler equivalent width measurement given by PAHFIT. For metal-rich systems, plot- 
ting these two ratios against one another results in a clear distinction between star-forming 
galaxies (stars and squares) and AGN (triangles), the latter having a higher [Olv]/[Nell] 
ratio and a weaker 7.7 /tm PAH feature than the former. However, the situation for the 
BCDs in our sample is different. They do not fall in the star-forming but rather cl o ser to 
the AGN region. The starburst with the highest [O iv] fraction in the IGenzel et al.l (119981 ) 
sample is NGC525 3 (open star symbol), a BCD-like galaxy with an oxygen abundance of 
12+log(0/H)~8.2 ( IKobulnicky et al.l Il999l ). Two BCDs in our sample would have been 
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classified as pure AGN on the basis of this diagnostic diagram. Thus, the strength of [O iv] 
combined with the 7.7 /xm PAH deficit in metal-poor starbursts makes them difficult to 
distinguish from AGN. 



The sec ond AGN diagnostic d iagram often used is illustrated in Fig. |20j This diagnostic, 
proposed by lLaurent et al.l (120001 ). exploits the steeply rising MIR continuum in starbursts, 
compared to the flatter continua in AGN and quiescent star- forming regions (or PDRs). 
Most star-forming galaxies have both PDR and Hll region components in their spectra; the 
PDR component tends to have a flat MIR continuum and be dominated by PAH bands, 
while the Hll regions have a deficiency of PAHs and a steeply rising MIR continuum from 
warm dust. To analyze our BCD sam ple, we have not used the original diagnostic diagram 
as proposed by Laurent et al.l (120001 ). but a slightly modified version of it as formulated 
by lPeeters et al.l (j2004al ). It differs from the original diagram in the continuum wavelength 
range that enters in the denominator of the labels on both axes. Fig. [20] plots the 14—15 /im 
continuum against the 6.2 /im PAH feature, both normalized by the 5.3—5.8 /im continuum. 
In this diagram, because of their steep MIR continua, the BCDs (filled circles) are clearly 
distinguished from the AGN (triangles), ULIRGs (squares), and metal-rich star- forming 
galaxies (stars), lying above them. A few BCDs are consistent with metal-rich Galactic Hll 
regions (x), but others are more extreme, both in their steeper continuum slope and in 
their enhanced PAH strength (or faint 5 /im continuum). The four BCDs with the steepest 
continua and most enhanced PAH strength are also the most metal-rich galaxies in our 
sample (Haro3, Mrk450, Mrkl 315, and UM311). W e conclude that the AGN/starburst 
diagnostic diagram proposed by lLaurent et al.l (120001 ) based on the continuum MIR slopes 
of these types of objects effectively separates AGN from starbursts, even at low metallicity. 



Molecules 



A long-standing puzzle of metal-poor star formation has been the apparent lack of cool 
molecular gas from which to form stars. Numerous attempts to detect carbon monoxide in 



< 



vigorous on- 


2;oing star formation, there seems to be very little CO in BCDs (t 


Saee et al. 


1992; 


Tavlor et al. 


1998; 


Gondhalekar et al. 


1998: 


Barone et al. 


2000: 


Lerov et al. 


2005 


). Intrigu- 



ingly, 12+log(0/H) £ 8.2 appears to be the same metallicity "threshold" below which PAH 



emission is thought to be supp ressed (jEngelbracht et al.ll2005l ; IWu et al.ll2006l ; iMadden et al. 



20061 ; lEngelbracht et al.l 120081 ) . Here we show that warm molecular gas in the form of H2 



does exist at very low metallicity. Our new IRS spectra reveal a variety of H2 rotational 
lines, and more than a third of the objects in our sample (8 BCDs) have ^ 3a detections in 
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one or more of the four lowest-order transitions of H2. 



6.1. Ho Emission 



Rotational transitions of H2 are important diagnostics of the warm neutral phase of 
the ISM, at temperatures between ~100 and 1000 K. These transitions, observable from 
space between 5 and 28 /xm, are a main coolant of the warm molecu lar component. Massive 
stars are the main source of excitation ( IHollenbach &: Tielenslll997l ). as H2 molecules can be 
pumped by non-ionizing far-ultraviolet (FUV) photons into excited states. They then return 
to the ground state through fluorescence or collisional deexcitation. The pure fluorescence 
mechanism is likely to be responsible for the near-infrared roto- vibrational transitions (e.g., 
Puxley et al.lll988[ ). because of the higher critical densities needed for collisional deexcitation. 
Critical densities of the lower-order rotational transitions are sufficiently low ( £ 10 3 cm -3 for 
S(0) to S(3)) that the lower rotational levels should be thermalized. For the pure rotational 
H 2 transitions studied here, Av is 0; we will refer to the AJ series of rot ational transitions 
as S(0), S(l), S(7) (for a complete tabulation, see lRoussel et al.ll2007l ). 



We first examine t he H2 emission i n term s of its contribution to the energy budget in the 
neutral ISM. Following IRoussel et al.l (120071 ). we use the sum of the S(0) to S(2) transitions 
to assess H 2 cooling. The bulk of H 2 cooling occurs through these lines. Figure [5H shows 
the sum of the S(0) to S(2) H 2 emission, normalized to the total IR luminosity, Lttr. As 
before, both the SINGS and BCD samples are plotted; the SINGS galaxies are confined to 
the hatched region. The S [Ho, (0 — 2 )1 / TIR ratios range from 1 x 10~ 4 to 2 x 10 -3 , with a mean 
of ~ 5 x 10~ 4 for the BCD sa mple R The E[H 9 f0 - 2)1/TIR for the SINGS galaxies range 
from 2.5 x 10~ 4 to 7.5 x 10~ 4 (IRoussel et al.l 120071 ) . Relative to Ltir, the warm molecular 



content of the BCDs is comparable to that of the SINGS galaxies. 

Because warm H 2 emission originates in PDRs, the UV-illuminated parts of molecu- 
lar clouds, H 2 and PAH emission should be related, since both species are excited by UV 
photons. Although our sample contains too few points to explore such a correlation, the 
mean ratio of S[ H 2 (0 — 2)| to S(PA H) of ~6% is roughly ten times larger than that for the 
SINGS galaxies (IRoussel et al.l 120071 ). Because the H 2 emission of BCDs, relative to TIR, is 
similar to the SINGS sample, this difference arises mainly from the deficit in PAH emission 
at low metallicites. Models show that the i ntensity of rotational transit ions of H 2 peaks 
inside molecular clouds at extinctions Ay~2 (IHollenbach fc Tielenslll997l ). while PAHs are 
expected to be excited predominantly on the cloud surface layers. Deeper into the cloud, 



'We have excluded UM311 because of the probable overestimate of its IR flux. 
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densities are higher and temp eratures lower, and PAHs would be coagulated onto grain 
mantles (jBoulanger et al.lll990f ). The similarity of H 2 /Ltir over a wide range of metallicity 
implies that H 2 might be somewhat self-shielded in BCDs. Although PAHs can survive the 
impact of FUV photons with energies between ~11 and 13.6 eV, H 2 would be dissociated 
without self-shielding. It is also difficult to understand how, without self-shielding, the H 2 
could survive the intense UV radiation field in low metallicity BCDs that is destroying the 
smallest PAHs (see floUj) . 

Although more than a third of our sample show definite low-order H 2 emission in their 
IRS s pectra, all attempt s to detect H 2 in BCDs in UV FUSE absorption spectra0 have 
failed ( IThuan et al.ll2005l ). The absence of H 2 absorption lines in FUSE spectra implies that 
the warm H 2 detected through IR emission must be quite clumpy. The FUSE observations 
are not sensitive to such a clumpy H 2 distribution because they can only probe the trans- 
parent UV sight lin es, not being able to penetrate dense clouds with N(H 2 ) ~ 
faoopes et alll2004h . 



10 20 cm- 2 



6.1.1. H 2 Excitation Temperatures 

Excitation diagrams for the H 2 lines help infer the excitation temperatures for molecular 
hydrogen. For four galaxies in our sample, there were sufficient data to construct excitation 
diagrams. These diagrams show, as a function of the upper level energy fcTj, the natural 
logarithm of the column density of the species N in the upper level of each transition, divided 
by the statistical weight, Assuming local thermodynamic equilibrium (LTE), the total 
column density N tot then follows from the Boltzmann equation: 

exp -— (1) 



Atot Z(T ex ) \ T ex 



where T ex is the excitation temperature for the i th level, and Z(T ex ) is the partition function 
The statistical weight is given by gi = {21 + 1)/ (2 J + 1); the spin number ta kes the values 



1 = for even J (para) transitions, and 1 = 1 for odd J (ortho) transitions. iHerbst et al. 



(119961 ) gives a convenient expression for the partition function: Z(T ex ) = 0.0247 T/[l — 
exp(— 6000 K/T)], valid for temperatures ^ 40 K. The column density for a given transition 
Ni is calculated from the observed H 2 flux Fi in that transition: 

FiAix 

1 hvjAM ^ 



Half a dozen have been observed, none of which are in the present sample. 
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where A{ is the transition's spontaneous emission coefficient, i>i is the frequency of the tran- 
sition, and Q is the beam solid angle. 

Figure I2"2l shows the resulting excitation diagrams for the four BCDs which have at least 
one detection at the 3cr level or greater in a low-order transition (except for SBS 1152+579 
for which we have used 2.5a detections), and at least two other significant ( ^ 3a) detections 
in higher-order lines. Since the inferred column density for a given transition depends on 
the beam size, we have experimented using different beam sizes of the IRS modules for the 
various transitions. The main result of this exercise is that the most appropriate beam size 
is the largest LH beam, equal to ll / .'lx22 / .'3. Because the scale of our spectra was ultimately 
set by the LH background subtraction (see SjSJ), and it was never necessary to multiply by a 
scaling factor larger than the ratios of the apertures, this is a reasonable conclusion. However, 
the sources observed at 24 jum with MIPS are either unresolved or only slightly extended, 
which would imply a slightly different beam dilution for each galaxy at longer wavelengths; 
applying such a correction would have been highly problematic. Nevertheless, since we have 
used the IRS point-source flux calibration, this may not be a serious problem, and, because 
of the way our spectra have been scaled, the excitation diagrams should not be affected 
significantly by this approximation. With a beam size that is too small, column densities 
are overestimated and total masses are underestimated. Hence, in Fig. |22J we have taken 
the most conservative approach for column densities, using the LH beam as discussed above. 
At the median distance of 21Mpc for our sources, this corresponds to a region of ~1.6kpc 
in diameter. 

The excitation temperature of the line-emitting gas is the reciprocal of the slope of the 
excitation curve; in the absence of fluorescence, this temperature would be the gas kinetic 
temperature in LTE. Assuming thermal emission^!. Fig. [22] implies a range of temperatures 
for the molecular gas. As a crude approximation to the probable continuous temperature 
distribution of the gas (see the parabolic fit for CGC G005 — 027), we have fitted the data 



points in each diagram by two linear segments (see also lRigopoulou et al.ll2002t iHigdon et al. 



20061 ; iRoussel et al.l 120071 ). The temperature inferred is a strong function of the set of tran- 
sitions from which it is derived: lower-order transitions probe lower-temperature gas than 
the higher-order ones. The implied temperature for the warm gas (from the lower-order 
transitions) is ~245K (not considering Mrk996 at 98 K, the only galaxy with a significant 
S(0) detection). The temperatures of the hotter molecular component (from the high-order 
transitions) range from ~820 K to ~1600 K. Considering the transitions from which they are 



12 This assumption is almost certainly invalid for the higher-order transitions, because that would require 
densities XL 10 5 cm -3 . In any case, the excitation diagrams remain a useful tool for inferring physical 
conditions of the line-emitting gas. 
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derived, these temperat ures are consistent w ith those found for starbursts (jRigopoulou et al. 
2002h . SINGS galaxies (IRoussel et al.ll2007h . and ULIRGs Jrligdon et alJbood l. 



We have not considered any departure fr om the standard ortho-para ratio of 3 for LTE 



(e.g., iBurton et al.lll992l ; iRoussel et al.ll2007l ). Departure from thermalization of ortho and 
para levels would be suggested by a non-monotonic increase of inferred excitation temper- 
atures with increasing rotational transition. The temperatures derived from mixed-parity 
H 2 transitions for the BCDs in our sample are generally consistent, indicating no departure 
from thermalization for most objects, although Fig. [22] does show some sign of a discrepancy 
for SBS 1152+579. A detailed discussion of this topic at low metallicity is precluded by the 
paucity of the BCD data presented here. 



6.1.2. Column Densities and Gas Masses 

Once we have an estimate of the gas temperature for each galaxy, we can derive the H 2 
column density and total H2 mass within the IRS beam. The derived masses and column 
densities depend on the size of the IRS beam used to relate the observed flux to the column 
densities. As stated above, we adopted a conservative approach for the column densities, 
and used the large LH beam, which may give masses that are slightly overestimated. 

In our discussion of H2 column densities and masses, we will consider only those objects 
in which "secure" temperatures have been derived, i.e. those that have two relatively reli- 
able detections (at the 2.5<r level or greater) of lower-order transition^. Four objects meet 
these criteria: CGCG005 - 027, HS 0837+4717, Mrk996, and SBS 1152+579. Mrk996 and 
CGCG005 - 027 have a total warm (T~100-120K) H 2 column density N H2 ~3x 10 21 cm" 2 
or 53M pc -2 . At a higher temperature of ~340K, SBS 1152+579 has a much lower molec- 
ular col umn density with Nh2 7 ^ 3 x 10 18 (0.04 M , pc~ 2 ). With a warm/total H 2 fraction 



of ~5% (jRigopoulou et al.ll2002t IRoussel et al.ll2007l ). the total molecular column density in 



SBS 1152+579 would be £ 1 M pc~ 2 . The column density for HS 0837+4717 is more uncer- 
tain, because it has only 2.5cr detections; the inferred N H2 is high, ~ 5 x 10 21 (78M pc~ 2 ), 
at a temperature of ~95K. With the exception of SBS 1152+579, the inferred H 2 column 
densities derived for the BCDs are rather high, compared to the values normally thought to 
hold for metal-poor galaxies. 

Converting to masses, the two most extreme cases are Mrk 996 with a total warm (98 K) 
H 2 mass in the IRS beam of 1.5xlO 8 M and CGCG005 - 027 (at 120 K) of 1.3xlO 8 M . 



13 Now including S(3) to consider correctly parity in the derivation of temperature. 
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Compared to the SINGS galaxies, these are rather high warm H2 masses. There are only 
three SINGS galaxies with warm ( £ 100 K) H 2 masses which exceed 10 8 M Q , and only one 
(NGC7552) with a warm H 2 mass greater than that of Mrk996. The other BCDs in our 
sample, with at least one low-order detection have warm (250-470 K) H 2 masses ranging from 
4 x 1O 3 M (Mrk209) to 5 x 1O 7 M (HS 0837+4717). T he median SINGS w arm H 2 mass 
is ~ 3.4 x 1O 6 M , and ranges from ~ 10 3 to 3 x 1O 8 M jRoussel et alil2007h . Most of the 
BCDs have £ 10 3 times more warm H 2 m ass than the SINGS low-metallicity galaxies with 
CO detections, NGC6822 and NGC2915 faoussel et aDl2007h . 



As stated above, these mass estimates may be overestimates. The true beam size is 
probably smaller than that of LH. Moreover, molecular hydrogen is not distributed uniformly 
within that beam but in dense clumps. However, even if we reduced the mass estimates by a 
factor of five (roughly the difference in beam area between LH and SH), the BCDs' molecular 
content would still be quite high. This is contrary to the conventional view that H 2 molecules 
do not exist in a low-metallicity environment, because of the lack of grains on which to form 
them. Clearly, metal abundance is not the only factor guiding H 2 formation in low-metallicity 
BCDs. 

Finally, we can compare the warm H 2 mass to that of the neutral atomic gas. According 
to the star formation models by lKrumholz et al.l (120091 ). the H 2 column density of Mrk966 
would correspond to an ISM m olecular fraction ~5 0%, given its metal abundance. Mrk996 
has an Hi mass of 1.5x 10 8 M (jThuan et al.lll999l ); thus its H 2 warm mass is comparable 
to that of its n e utral atomic gas mass, in agreement with the theoretical predictions of 
Krumholz etaD J2009h . 



Surprisingly, our data suggest that metallicity and H 2 content are not directly linked. 
Of the four BCDs discussed above with significant low-order H 2 detections, three are below 
the sample median 12+log(0/H) of 7.9. Other factors must play a role in determining the 
H 2 fraction. In star-forming regions where substantial amounts of dust are concentrated in 
high-density regions, self-shielding can promote H 2 formation; thus the re gion can cool more 
effectively and star- formation rate is enhanced (IHirashita &: Huntl 120041 ). Thus H 2 content 
may be more correlated with compactness (size and density) than with metal abundance. 



6.2. Water and the Hydroxyl Radical 

Our IRS spectra tentatively suggest the presence of H 2 and OH in six objects (see 
Table [6]). Figure [3] shows the ^ 3cr detections of H 2 in three BCDs and OH in four. The 
H 2 lines at 29.836 /zm and 29.885 fim result from rotational transitions of H 2 16 7 2 5 — 6ig 
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and 542 — ^13, respectively. The OH emission lines correspond to rotationally excited levels 
2 n 1/2 -> 2 n 3/2 7/2 - 5/2 (28.940 /xm), 2 n 3/2 19/2 - 17/2 (30.346 /xm), and 2 n 1/2 17/2 - 15/2 
(30.657 /im). Because some of the spectra also show spurious spectral features near the 
tentative detections, not associated with known emission lines, these molecule detections are 
uncertain. Nevertheless, if real, the transitions of "hot" H 2 and OH would be, to the best 
of our knowledge, the first such detections in extragalactic objects. 

Molecular emission from H 2 and OH is perhaps not unexpected, given the physical 
conditions in some BCDs. H 2 e mission from vapor is thought to arise from slow nondis- 
sociative, or C-type shocks (e.g., iDraind Il980l ). In such a situation, at a shock velocity 



of £ lOkms -1 , temperatures would exceed 300 K, and s ubsequent reactions cou ld rapidly 
convert all gas-phase oxygen not bound in CO into H 2 (jElitzur &: de Jonglll978l ). If there 
is a connection between the dissociative shocks responsible for H 2 vapor emission and the 
fast radiative shocks possibly responsible for [Oiv] emission (see § 14.41) . we would expect 
both types of emission. Indeed, all three galaxies with H 2 emission in our sample also 
show £ 3er [O iv] detections. Nevertheless, to be consistent with the shoc k scenario for H 2 Q 
prod uction, the H 2 density would need to be rather high, £ 10 5 cm -3 (jElitzur fc de Jong 
19781 ). The electron densities in the ionized gas inferred from the [S ill] line ratio (§ I4.3P are 
considerably smaller than this, but the [Oiv] and H 2 emission could arise from different 
regions that are much denser. We have no way of verifying this condition with the present 
data. With Herschel, it may be possible to further pursue observations of H 2 vapor at low 
metallicity, and better understand the necessary physical conditions and formation scenarios 
for this molecule. 

OH emission could be associated with the photo- dissociation of H 2 Q at photon energies 
£ 9eV, similar to the situation in Galactic outflows ( iTappe et al.l 120081 ) . Interestingly, the 
four galaxies with OH detections also have high-order H 2 detections, implying the presence 
of hot dense ga s, thought to be necessary f or the neutral reactions leading to OH and H 2 



formation (e.g.. iHollenbach &: McKedll979i ) 



7. Summary and Conclusions 

We have presented low- and high-resolution IRS spectra, supplemented by IRAC and 
MIPS measurements, of 22 BCD galaxies, obtained during our Cycle 1 GO Spitzer program 
(PID 3139). The BCD sample was chosen to span a wide range in oxygen abundance 
[12+log(0/H) between 7.4 and 8.3], and in ISRF hardness as measured by the intensity of 
the nebular Hell A 4686 emission line relative to H/3. The IRS spectra provide a variety of fine- 
structure lines, aromatic features, and molecular lines which probe the physical conditions 
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in a metal-poor ISM, and enable a study of the dust properties as function of metallicity, 
hardness and intensity of the ISRF. We have used the PAHFIT routine to fit simultaneously 
the spectral features, the underlying continuum and the extinction in the IRS spectra. Fluxes 
have also been derived by fitting Gaussian profiles to all FS and molecular lines. To place the 
results for our low-metallicity BCD sa mple in perspective, w e have compared its properties 
to those of the SINGS galaxy sample (IKennicutt et al.ll2003l ). 

We have obtained the following results: 



(1) The [Ne in] /[Neil] and [Siv]/[Siii] flux ratios are good diagnostics for the softer UV 
radiation (< 40 eV). The [Ne lll]/[Ne il] ratio depends on metallicity, being 1-2 orders of 
magnitude greater in our metal-poor BCD sample than in the more metal-rich SINGS 
galaxies, but the ratio flattens out at 12+log(0/H) £ 8.3. The [Olv]/[Sill] ratio is an 
effective measure of harder radiation, and there is a strong correlation of the [O lv]/[Si il] 
ratio with oxygen abundance, implying that lower metallicity galaxies have harder ion- 
izing radiation than more metal-rich ones. 

(2) We detect at 3a [Oiv] (ionization potential of 54.9 eV) in 7/22 observed galaxies, but 
[Fell] (ionization potential of 7.9 eV) in only 2. At low metallicity, [Oiv] emission is 
almost 4 times as common as [Fell] (32% vs. 9%). This is another indication that the 
ISRF of low-metallicity galaxies is very hard. 

(3) Electron densities derived from the IR [S ill] line ratio do not correlate with those in- 
ferred from the optical [S il] ratio, or with the presence of [O iv] emission. We would 
expect a correlation of [O iv] intensity with electron density, if fast radiative shocks were 
responsible for the ionizing radiation that produces [O iv] . Perhaps the [S ill] IR lines 
arise in and probe a region that is less dense than the [O iv] region. 

(4) The ratio of 71 to 160/mi fluxes (z/Lv(71))/ (z/L^lGO)) is sensitive to the temperature of 
the large ("classical") grains, and is thus a good indicator of the ISRF intensity. The 
latter is ^ 5 times greater in the BCDs than in the SINGS galaxy with the most intense 
ISRF, implying that metal-poor star-forming galaxies have not only a harder ISRF, but 
also a more intense one than metal-rich galaxies. 

(5) Two-thirds of the BCDs show PAH features. The flux ratios of the different bands are 
typical of the largest PAHs modeled so far, with iV min £ 100 C atoms. We interpret these 
trends as an indication that low-metallicity BCDs contain relatively larger PAHs than 
more metal-rich environments. Apparently, only these large PAHs are able to survive 
the hard and intense ISRF in a low-metallicity ISM. 
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(6) The fraction S(PAH)/TIR of PAH emission S(PAH) normalized to the total IR lumi- 
nosity TIR is considerably smaller in low-metallicity BCDs (~0.5%) than in the more 
metal-rich SINGS galaxies (~10%). At low metallicity, just as in metal- rich galaxies, 
the PAH fraction of Ltir is constant. 

(7) There is a good correlation between S(PAH)/TIR and the [Ne ill]/ [Ne n] flux ratio, but 
not with metallicity or with ISRF intensity. Evidently, over the range in UV energies 

£j 41 eV, the fraction of TIR emerging as PAH emission depends on the hardness of 
the ISRF. This suggests that the PAH fraction in BCDs is not directly controlled by 
metallicity, but rather by the hardness of the radiation field, also responsible for the 
destruction of the smallest PAH particles. 

(8) The hardness of the ISRF in BCDs can be comparable to that in AGN. Because of 
this, the [Olv]/[Nell] flux ratio, often used as a diagnostic to distinguish star-forming 
galaxies from AGN at solar abundances, cannot play that role in the low-metallicity 
regime: in the [Oiv]/[Neii] vs. 7.7 /mi PAH band diagnostic diagram, low-metallicity 
BCDs occupy the same region as metal-rich AGN. On the contrary, the AGN/starburst 



diagnostic diagram proposed by ([Laurent et al.l 120001 ) . based on the continuum MIR 



slopes of these types of objects, does a good job at separating AGN from starbursts, 
even at low metallicity. 

(9) Our IRS spectra reveal a variety of H 2 rotational lines, and more than one third of 
the objects in our sample (8 BCDs) have ^ 3a detections in one or more of the four 
lowest-order transitions of H2. The BCDs contain more H2 than most SINGS galaxies, 
relative to Ltir- This is quite contrary to the usual assertion that H 2 molecules do 
not exist in a low-metallicity environment, because of the lack of grains on which to 
form them. Clearly, metal abundance cannot be the only factor guiding H 2 formation in 
low-metallicity BCDs. 

(10) The mean ratio of H 2 to PAH emission is ~6%, or 10 times larger than that for the 
SINGS galaxies (but comparable to SINGS AGN). This difference arises mainly from 
the BCD deficit in PAH emission. While PAHs can survive the impact of FUV photons 
with energies between ~11 and 13.6 eV, H 2 would be dissociated without self-shielding, 
implying that H 2 could be somewhat self-shielded in BCDs. 

(11) The mean excitation temperature derived for the warm gas from the lower-order H 2 tran- 
sitions is ~245 K. The temperatures of the hotter molecular component from the high- 
order H 2 transitions range from ~820K (in Mrk996 and SBS 1152+579) to ~1600K. 
These temperatures are consistent with those found for starburst SINGS galaxies and 
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ULIRGs. The warm molecular eras masses in our BCDs range from 10 3 to 1O 8 M , and 
can be comparable to the neutral hydrogen gas mass. 

(12) Some of our IRS spectra suggest tentatively the presence of molecules other than H 2 , 
such as H2O and OH. All three galaxies with H2O emission also show a ^ 3a [O iv] line, 
consistent with a shock scenario for H 2 production. All four galaxies with OH emission 
also have high-order H 2 detections, implying the presence of hot, dense gas which would 
favor OH and H 2 formation. 



Our IRS data have shown that the ISRF in metal-poor BCDs can be as extreme in 
hardness as the radiation in AGN. In fact, the deficit of PAH emission at low metallicity and 
their hard radiation as indicated by [O iv] lumin osity would have plac ed some BCDs as AGN 
in diagnostic diagrams such as that proposed by lGenzel et al.l (119981 ). Our analysis suggests 
that in the absence of an AGN, low metal abundance is a necessary condition for such 
extreme fields. However, we have also shown that low metallicity alone does not guarantee 
an extreme ISRF in terms of hardness and intensity. Other factors such as compactness of 
the star-forming region must play a role. 

We have inferred a deficit of small PAHs, and find that PAH destruction mechanisms 
are more closely related to the hardness of the radiation field, rather than to metallicity. 
This is a rather difficult exercise, since most of the [O iv] detections are in galaxies where we 
have little information on the PAH population, and the converse is also true. Nevertheless, 
there appears to be a close connection between hardness of the ISRF and the properties of 
the PAHs; this connection is tighter than that with metallicity. In fact, we find correlations 
of PAH properties with ISRF hardness, but almost none with metallicity, in spite of the 
good correlations between hardness tracers and metal abundance. 

It is not clear whether the origin of such hard radiation can be stellar. In their de- 
tailed modeling of the [O iv | line in Mrk 996, iThuan et al.l ( 120081 ) ran Costar models by 
Schaerer fc de Koter I (119971 ) using stars with the highest effective temperature (T e ff =53,000 
K) and the hardest radiation possible (corresponding to the hottest 03 stars). They failed 
to reproduce the observed intensity of the [O iv] line by a factor of seve ral. They then 



consid ered Wolf-Rayet stars of type WNE-w. According to calculations by ICrowther et al. 



(119991 ). models for WNE-w stars (early nitrogen Wolf-Rayet stars with weak lines) show 
a strong ionizing flux above 54.4 eV, in contrast to the WCE and WNL stars, which show 
negligible fluxes above that energy. The line intensi ty of [O IV"1 in Mrk 996 can be reproduced 
with such a population of WNE-w stars. However, IThuan et al.l (120081 ) do not favor such an 
hypoth esis because WNE-w st ars are very rare. Rather, those authors favor fast radiative 
shocks (IThuan fc Izotovl 120051 ) as the origin of the high-ionization radiation. Shocks are a 
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more realistic expectation in dense compact regions with intense star formation because of 
SN explosions and massive stellar winds. At low metallicity, the less efficient cooling may 
intensify the effect of the shock. We postpone an in-depth discussion of the origin of the 
hard ionizing radiation for the whole of our sample to a future paper. There, we will present 
detailed photoionization and shock models to account for both optical+ IR emission line 
intensities. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by JPL/Caltech under NASA contract 1407. Support for this work was provided by 
NASA Spitzer GO grant JPL-1263707. We acknowledge financial contribution from contract 
ASI-INAF 1/016/07/0. 
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Table 1. Sample Galaxies 



Name 


Distance a 


Redshift 


12+log(0/H) 


c H p h 


Heii/H/3 


CGCG007-025 


20.8 


0.00483 


7.76 


0.29 


0.0136 


HS0837+4717 


174.3 


0.04195 


7.60 


0.29 


0.0228 


HS1442+4250 


12.6 


0.00211 


7.63 


0.11 


0.0297 


HS2236+1344 


86.4 


0.02062 


7.47 


0.16 


0.0106 


Haro3 


17.4 


0.00323 


8.32 


0.24 


0.0000 


IIZw70 


21.8 


0.00394 


8.04 


0.21 


0.0056 


J0519+0007 


181.8 


0.04476 


7.43 


0.30 


0.0249 


Mrk5 


15.4 


0.00264 


8.04 


0.42 


0.0387 


Mrk36 


7.6 


0.00215 


7.81 


0.02 


0.0244 


Mrk209 


4.8 


0.00094 


7.81 


0.00 


0.0117 


Mrk450 


15.1 


0.00288 


8.15 


0.14 


0.0031 


Mrk724 


20.0 


0.00402 


8.03 


0.12 


0.0000 


Mrk996 


21.8 


0.00541 


8.10 


0.53 


0.0000 


Mrkl315 


13.1 


0.00282 


8.25 


0.16 


0.0000 


Mrkl329 


13.1 


0.00544 


8.25 


0.16 


0.0000 


SBS0917+527 


35.4 


0.00776 


7.90 


0.09 


0.0221 


SBS0946+558 


25.5 


0.00517 


8.04 


0.18 


0.0149 


SBS1030+583 


35.3 


0.00757 


7.83 


0.00 


0.0234 


SBS1152+579 


74.5 


0.01720 


7.85 


0.26 


0.0152 


SBS1415+437 


11.7 


0.00203 


7.61 


0.00 


0.0227 


Toll924-416 


38.4 


0.00945 


7.94 


0.26 


0.0190 


UM311 


23.0 


0.00559 


8.31 


0.15 


0.0000 



a Distance in Mpc, taken from NED. 
b Extinction coefficient: A H( g=2.5,c_f//3- 
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Table 2. Short- Wavelength Fine-Structure Lines 3, 



Name 


[Ar n] 


[Ar in] 


[Siv] 


[Ne n] 




6.985/um 


8.991/im 


10.511/zm 


12.814^m 




15.759 eV 


27.629 eV 


34.830 eV 


21.564eV 


CGCG 007-025 


••• (0.33) 


0.66 (0.18) 


5.98 (0.43) 


••• (1.29) 


Haro3 


3.44 (0.97) 


10.60 (1.67) 


26.80 (3.70) 


27.50 (3.59) 


HS 0837+4717 


••• (0.13) 


0.28 (0.09) 


3.07 (0.41) 


••• (1.23) 


HS 1442+4250 


••• (0.21) 


••• (0.21) 


••• (1.60) 


••• (1.60) 


HS 2236+1344 


••• (0.12) 


••• (0.12) 


2.50 (0.56) 


••• (1.68) 


IIZw70 


••• (0.72) 


••• (0.72) 


4.03 (0.26) 


1.34 (0.30) 


JUoiy-t-uuu ( 


• • • \\).ZZ ) 


• • • [VJ.ZZ ) 


( 1 1 7\ 

••• (1.17) 


( 1 1 7\ 

••• (1.17) 


Mrk5 


0.54 (0.08) 


0.52 (0.08) 


1.64 (0.13) 


1.08 (0.13) 


Mrk 36 


••• (0.27) 


••• (0.27) 


2.36 (0.14) 


••• (0.43) 


Mrk209 


1.26 (0.33) 


1.58 (0.26) 


8.63 (0.18) 


••• (0.45) 


Mrk 450 


••• (0.21) 


1.46 (0.13) 


4.45 (0.35) 


0.91 (0.29) 


Mrk 724 


••• (0.28) 


••• (0.28) 


3.92 (0.44) 


••• (1.31) 


Mrk 996 


••• (1.23) 


1.87 (0.40) 


3.90 (0.32) 


4.55 (0.32) 


Mrk 1315 


0.73 (0.05) 


1.36 (0.12) 


10.20 (0.48) 


••• (1.43) 


Mrk 1329 


••• (0.46) 


1.80 (0.15) 


12.30 (0.27) 


1.27 (0.27) 


SBS 0917+527 


••• (0.12) 


••• (0.12) 


1.46 (0.12) 


••• (0.36) 


SBS 0946+558 


••• (0.15) 


••• (0.15) 


••• (1.63) 


••• (1.63) 


SBS 1030+583 


••• (0.10) 


••• (0.10) 


1.93 (0.59) 


••• (1.77) 


SBS 1152+579 


••• (0.40) 


0.92 (0.13) 


6.32 (0.12) 


••• (0.36) 


SBS 1415+437 


••• (0.18) 


••• (0.18) 


1.02 (0.22) 


••• (0.49) 


Tol 1924-416 


••• (0.80) 


2.56 (0.27) 


11.40 (0.34) 


1.88 (0.34) 


UM311 


••• (0.29) 


1.71 (0.23) 


5.22 (0.69) 


••• (2.08) 



a Fluxes are measured by integrating Gaussian profiles as fit by PAHFIT, 
and are in units of 10~ 17 Wm~ 2 . la uncertainties are given in parentheses; 
when there is no flux available, the value reported in parentheses is the 3a 
upper limit. 



Table 3. Dust Features a 



Name 


DF_5 


DF_6 


DF_7 b 


DF.8 


DF.ll 


DF_12 


DF_17 


SPAH 




5.7/zm 


6.2/xm 


7.7/un 


8.6^m 


11.3^m 


12.6^m 


17^m 


Total c 
















Complex 




CGCG 007-025 


••• (0.33) 


••• (0.33) 


6.34 (0.82) 


••• (0.33) 


2.81 (0.51) 


••• (1.28) 


••• (0.33) 


12.79 (1.12) 








7.459 




11.199 














0.470 




0.067 








Haro3 


20.30 (3.53) 


101.00 ( 10.60) 


188.50 ( 31.85) 


56.30 ( 14.40) 


148.90 ( 41.95) 111.00 ( 19.30) 


••• (1.39) 680.60 ( 60.62) 




5.744 


6.245 


7.808 


8.682 


11.335 


12.673 








0.188 


0.194 


0.362 


0.279 


0.544 


0.798 






HS 0837+4717 


••• (2.91) 


••• (0.13) 


••• (0.13) 


••• (0.13) 


1.37 (0.41) 


••• (1 23) 


••• (0 13) 


2.95 (0.57) 












11.214 


















o non 

u.uyu 








HS 1442+4250 


••• (0.21) 


0.94 (0.24) 


••• (0.21) 


••• (0.57) 


3.46 (0.53) 


••• (1.60) 


••• (0.21) 


6.15 (0.79) 






6.294 






11.347 












0.094 






0.545 








HS 2236+1344 


••• (0.12) 


••• (0.12) 


••• (0.12) 


••• (0.12) 


••• (1.68) 


••• (1.68) 


••• (0.12) 


••• (0.12) 


IIZw70 


••• (0.72) 


••• (3.78) 


10.81 (1.98) 


••• (0.72) 


••• (0.66) 


••• (0.66) 


••• (0.72) 


12.39 (2.03) 








7.808 


















0.207 












J0519+0007 


••• (0.22) 


••• (0.22) 


2.45 (0.55) 


••• (0.22) 


■•• (1.17) 


••• (1.17) 


••• (0.22) 


4.24 (0.67) 








7.405 


















0.467 












Mrkl315 


••• (0.15) 


3.97 (0.12) 


1.91 (0.10) 


3.47 (0.48) 


1.45 (0.48) 


••• (1.43) 


••• (0.15) 


15.24 (0.75) 






6.201 


7.500 


8.648 


11.256 












0.279 


0.165 


0.431 


0.088 








Mrk 1329 


••• (5.70) 


••• (3.03) 


9.31 (0.99) 


4.70 (0.87) 


6.94 (0.54) 


1.89 (0.27) 


1.85 (0.49) 


24.69 (1.53) 








7.520 


8.537 


11.272 


12.627 







Table 3 — Continued 



Name 


DF.5 
5.7^m 


DF_6 
6.2/xm 


DF.7 b 
7.7/«m 


DF.8 
8.6^m 


DF.ll 
11.3^m 


DF_12 
12.6/im 


DF_17 
17^im 
Complex 


SPAH 
Total c 








0.474 


0.166 


0.541 


0.265 






Mrk209 


••• (2.82) 


••• (0.32) 


••• (2.37) 


••• (0.32) 


••• (0.45) 


••• (0.45) 


••• (0.32) 


7.41 (1.28) 


Mrk36 


••• (3.00) 


••• (0.27) 


••• (0.27) 


••• (0.27) 


••• (0.43) 


••• (0.43) 


••• (0.27) 


4.73 (1.32) 


Mrk450 


2.06 (0.25) 


6.16 (0.23) 


16.40 (1.42) 


5.95 (0.85) 


5.46 (1.61) 


••• (0.76) 


••• (0.21) 


36.03 (2.33) 




5.800 


6.291 


7.520 


8.520 


11.331 










0.169 


0.283 


0.659 


0.255 


0.388 








Mrk5 


2.17 (0.20) 


3.27 (0.14) 


7.47 (0.71) 


3.45 (0.39) 


5.42 (0.64) 


••• (0.39) 


••• (0.24) 


24.33 (1.08) 




5.791 


6.287 


7.801 


8.708 


11.326 










0.214 


0.283 


0.281 


0.440 


0.544 








Mrk 724 


••• (0.28) 


••• (0.28) 


••• (0.28) 


••• (2.16) 


••• (1.31) 


••• (1.31) 


••• (0.28) 


••• (0.28) 


Mrk996 


7.03 (1.45) 


15.50 (4.32) 


17.90 (4.35) 


••• (9.63) 


11.80 (2.02) 


3.53 (0.71) 


••• (0.66) 


60.66 (6.71) 




5.757 


6.236 


7.493 




11.276 


12.684 








0.101 


0.179 


0.472 




0.250 


0.082 






SBS 0917+527 


••• (0.12) 


••• (0.12) 


••• (0.12) 


1.42 (0.23) 


1.00 (0.29) 


0.37 (0.12) 


••• (0.12) 


2.79 (0.39) 










8.559 


11.230 


12.601 














0.167 


0.539 


0.082 






SBS 0946+558 


••• (0.15) 


••• (0.15) 


0.89 (0.28) 


••• (0.15) 


2.74 (0.54) 


••• (1.63) 


••• (0.15) 


12.11 (1.03) 








7.520 




11.230 














0.474 




0.539 








SBS 1030+583 


••• (0.10) 


••• (0.10) 


4.75 (0.40) 


0.67 (0.14) 


••• (1.77) 


••• (1.77) 


••• (0.10) 


5.42 (0.42) 



Table 3 — Continued 



Name 


DF_5 


DF.6 


DF_7 b 


DF_8 


DF.ll 


DF.12 


DF.17 


SPAH 




5.7/zm 


6.2/mi 


7.7^m 


8.6^m 


11.3^m 


12.6^m 


17^m 
Complex 


Total c 








7.520 


8.534 
















0.474 


0.166 










SBS 1152+579 3.30 (0.70) 


••• (0.40) 


12.20 (0.67) 


3.41 (0.46) 


••• (0.36) 


••• (0.36) 


••• (0.40) 


24.15 (1.11) 




5.800 




7.520 


8.710 












0.305 




0.474 


0.510 










SBS 1415+437 


••• (0.18) 


••• (0.18) 


1 O O A /A OC\ 

12.84 (0.25) 
7.782 
0.258 


/ o a a \ 

■■■ (2.04) 


/a a n\ 

• • • (0.49) 


/A A A\ 

• • • (0.49) 


/A 1 0\ 

••• (0.18) 


12.84 (0.25) 


Tol 1924-416 


••• (0.80) 


••• (0.80) 


6.39 (1.40) 
7.507 
0.473 


2.08 (0.62) 
8.510 
0.166 


2.32 (0.58) 
11.230 
0.180 


3.09 (0.41) 
12.643 
0.266 


••• (0.80) 


13.88 (1.69) 


UM311 


••• (2.37) 


12.70 (2.14) 
6.216 
0.163 


38.76 (5.37) 
7.822 
0.207 


4.27 (1.29) 
8.634 
0.196 


5.53 (0.69) 
11.237 
0.067 


••• (2.08) 


••• (0.29) 


64.98 (6.00) 



a Dust features are integrated over the entire best-fit Drude profile as performed by PAHFIT. In the first line, 
fluxes are in units of 10~ 17 Wm~ 2 . la uncertainties are given in parentheses; when there is no flux available, the 
value reported in parentheses is the 3cr upper limit. The PAHFIT best-fit central wavelengths in ^im are given in 
the second line, and the best-fit FWHM in [im in the third. 

b This complex includes the sum of features from 7.4 to 7.8 ^m; the FWHM corresponds to the longest wavelength 
feature. 

c This total includes the sum of features in the table, together with features identified as PAHs by PAHFIT with 
wavelengths from 6.6 to 6.8 /xm, 8.3 ^m, and 10.6 to 10.8 /xm. 



Table 4. Long- Wavelength Fine-Structure Lines 3, 



Name 


[Ne in] 


[S nil 

[kj XXXj 


[Oiv] 


[Fen] 


[Fe in] 


[S ml 


ISi ill 


[Ne in] 




15. 555/i m 


18.713/um 


25.890/x.m 


25.988/um 


33.038Atm 


33.481^m 


34.815^m 


36.014/im 




40.962 eV 


23.330 eV 


54.934 eV 


7.870 eV 


16.180cV 


23.330 eV 


8.151 eV 


40.962 eV 


CGCG 007-025 


5.55 (0.49) 


••• (1.28) 




(1.49) 




(1.49) 


••• (1.77) 


3.48 (0.67) 


2.89 (0.50) 


••• (1.49) 


Haro3 


79.58 (0.27) 


40.99 (0.27) 


5.38 


(0.93) 




(2.80) 27.13 (0.93) 


116.41 (0.93) 


53.47 (0.93) 


12.11 (0.93) 


HS 0837+4717 


1.54 (0.41) 


••• (1.23) 




(0.92) 






... (0 Q9 , 


••• (0.92) 


••• (0.92) 


••• (0.92) 


HS 1442+4250 


••• (1.60) 


••• (1.60) 




(1.00) 




1 ± .UU I 


... (\ (\{\\ 

11 .UU J 


••• (1.00) 


••• (1.00) 


••• (1.00) 


HS 2236+1344 


••• (1.86) 


••• (1.68) 




(0.54) 




(ft 54i 


... (ft Z.A) 


••• (0.54) 


••• (0.54) 


••• (0.54) 


II Zw70 


6.78 (0.35) 


3.95 (0.22) 


0.77 


(0.17) 








10.09 (0.17) 


6.65 (0.53) 


••• (0.51) 


J0519+0007 


••• (1.17) 


••• (1.17) 




(0.67) 




(0.67) 


••• (0.67) 


••• (0.67) 


••• (0.99) 


••• (0.67) 


Mrk5 


4.31 (0.13) 


1.90 (0.13) 




(0.21) 


0.24 


(0.07) 


••• (0.21) 


5.02 (0.07) 


2.45 (0.07) 


1.01 (0.07) 


Mrk 36 


2.70 (0.14) 


1.53 (0.14) 


1.11 


(0.30) 




(0.91) 


••• (0.91) 


2.63 (0.30) 


2.97 (0.30) 


••• (0.91) 


Mrk 209 


5.57 (0.15) 


2.69 (0.15) 


1.46 


(0.19) 




(0.57) 


••• (0.57) 


5.66 (0.19) 


••• (0.57) 


••• (0.57) 


Mrk 450 


7.00 (1.33) 


4.27 (0.25) 




(0.64) 




(0.64) 


••• (0.64) 


5.68 (0.21) 


3.03 (0.21) 


1.96 (0.21) 


Mrk 724 


3.98 (0.44) 


2.89 (0.44) 




(1.73) 




(1.73) 


••• (1.73) 


••• (1.73) 


••• (1.73) 


••• (1.73) 


Mrk 996 


8.95 (0.32) 


7.31 (0.32) 


1.46 


(0.21) 


0.96 


(0.32) 


0.90 (0.21) 


16.43 (0.21) 


5.75 (0.21) 


••• (0.63) 


Mrk 1315 


9.74 (0.48) 


5.09 (0.48) 




(0.80) 




(0.80) 


••• (0.80) 


7.26 (0.27) 


1.81 (0.27) 


••• (0.80) 


Mrk 1329 


15.51 (0.27) 


8.07 (0.32) 




(2.10) 




(2.10) 


••• (2.10) 


8.59 (0.70) 


2.83 (0.70) 


••• (2.10) 


SBS 0917+527 


1.00 (0.12) 


0.67 (0.12) 




(0.76) 




(0.76) 


••• (0.76) 


••• (0.76) 


••• (0.76) 


••• (0.76) 


SBS 0946+558 


••• (1.63) 


••• (1.63) 




(0.51) 




(0.51) 


••• (0.51) 


1.56 (0.17) 


0.57 (0.17) 


••• (0.51) 


SBS 1030+583 


••• (1.77) 


••• (1.77) 




(0.83) 




(0.83) 


••• (0.83) 


••• (0.83) 


••• (0.83) 


••• (0.83) 


SBS 1152+579 


2.45 (0.12) 


1.38 (0.14) 




(0.76) 




(0.76) 


••• (0.76) 


2.68 (0.25) 


••• (0.76) 


••• (0.76) 


SBS 1415+437 


1.22 (0.16) 


1.09 (0.16) 


0.61 


(0.13) 




(0.40) 


••• (0.40) 


1.41 (0.13) 


2.15 (0.13) 


••• (0.40) 


Tol 1924-416 


23.99 (0.34) 


7.40 (0.34) 


3.33 


(0.61) 




(3.51) 


••• (1.84) 


18.57 (0.61) 


15.23 (0.74) 


••• (6.54) 


UM311 


13.75 (0.69) 


2.53 (0.69) 




(1.99) 




(1.99) 


••• (1.99) 


••• (1.99) 


4.01 (0.66) 


••• (1.99) 



a Fluxes are measured by integrating Gaussian profiles, and are in units of 10 17 Wm 2 . la uncertainties are given in 
parentheses; when there is no flux available, the value reported in parentheses is the 3cr upper limit. 
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Table 5. H 2 Lines a 
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(0.16) 
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(0.08) 




(0.76) 




(0.76) • • • 


(0.64) 




(0.66) 


1.15 


(0.35) 




(0.97) 


0.74 


(0.32) 0.71 


(0.21) 




(0.15) 




(0.15) 




(1.43) 




(1.43) ••• 


(0.80) 
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(0.18) 




(0.46) 




(0.80) 




(0.80) • • • 


(2.10) 


38 


(0.16) 




(0.12) 




(0.36) 




(0.36) • • • 


(0.76) 




(0.15) 




(0.15) 




(1.63) 




(1.63) ••• 


(0.51) 




(0.10) 




(0.10) 




(1.77) 




(1.77) ••• 


(0.83) 


41 


(0.15) 


0.39 


(0.13) 


0.30 


(0.12) 




(0.36) ••• 


(0.76) 




(0.18) 




(0.18) 


0.43 


(0.20) 




(0.49) • • • 


(0.40) 




(0.80) 


2.35 


(0.27) 




(1.03) 




(1.03) ••• 


(1.84) 


73 


(0.35) 




(0.29) 




(2.08) 


1.49 


(0.69) • • • 


(1.99) 



a Fluxes are measured by integrating Gaussian profiles, and are in units of 10 _17 Wm~ 2 . Unlike other tables, 
lines shown here as detections are 2cr. la uncertainties are given in parentheses; when there is no flux available, 
the value reported in parentheses is the 3a upper limit. 

b Unlike the shorter-wavelength lines which are fitted with Gaussian profiles by PAHFIT, these are measured by 



independent Gaussian profiles as described in the text. 
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Table 6. Water and Hydroxyl Lines a 



Name 


OH1/2-3/2 


o-H20 


P-H20 


OH19/2-17/2 


OH19/2-17/2 




28.939/um 


29.836^m 


29.885/zm 


30.346^m 


30.657,um 


Haro 3 


••• (2.80) 


••• (2.80) 


4.57 (0.93) 


••• (2.80) 


••• (2.80) 


IIZw70 


••■ (0.51) 


••• (0.51) 


••■ (0.51) 


0.59 (0.17) 


••• (0.51) 


Mrk5 


••• (0.21) 


••• (0.21) 


••• (0.21) 


••• (0.21) 


0.32 (0.07) 


Mrk996 


••• (0.63) 


2.48 (0.40) 


••• (0.63) 


••• (0.63) 


••• (0.63) 


Mrkl315 


1.65 (0.27) 


••• (0.80) 


••• (0.81) 


••• (0.80) 


••• (0.81) 


Tol 1924-416 


2.67 (0.61) 


••• (1.84) 


2.30 (0.61) 


••• (1.84) 


••• (1.84) 



a Fluxes are measured by integrating Gaussian profiles, and are in units of 10 -17 Wm~ 2 . 
la uncertainties are given in parentheses; when there is no flux available, the value reported 
in parentheses is the 3a upper limit. 
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Fig. 1. — IRS SL+SH spectra from 1 to 15/im, with the best-fitting PAHFIT model super- 
imposed shown as a red curve. The vertical axis is in units of mJy. Significant emission-line 
identifications are shown as labeled vertical lines. The top panels report the best-fit reduced 
xt, obtained over the wavelength region shown. The bottom panels show the residuals from 
the PAHFIT models. 



-45 - 




Fig. 1. — continued. 
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Fig. 1. — continued. 
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Fig. 2. — IRS LH spectra of the 26/xm region around [O iv] and [Fe il] for those objects with 
either [O iv] or [Fe il] detections. As in Fig. [TJ all marked features are considered significant 
(see Table H]). The MIPS24 total flux is marked with a red open circle. Only those spectral 
points with S/N>3 are shown. 
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Fig. 3. — Close-up of tentative H2O (29.8/im) and OH (28.9 jum) detections. All marked 
lines are considered to be significant detections (see Table [6]). Only those spectral points 
with S/N>3 are shown. 
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Fig. 4. — [Siv]/[Sm] (18.7 /im) flux ratios (top panel) and [Ne ill] /[Nell] ratios (bottom) 
plotted against the nebular oxygen abundance, 12 +log(0/H). The BCDs in our sample 
are shown as filled circles, while the SINGS galaxies (IDale et al.ll2009l ) as open squares (Hn 
region- nuclei) and open triangles (AGN). The horizontal dotted line in the lower panel shows 
the mean neon ratio (5.42) for 12+log(0/H)<8.3 (not considering lower limits). 
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Fig. 5. — [Siv]/[Sm] (18.7 /wn) flux ratios (top panel) and [Ne ill] /[Nell] ratios (bottom) 
plotted against the optical Hell/H/3 ratio. 
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Fig. 6. — [Oiv]/[Sin] flux ratios plotted against the optical Hell/H/3 ratio (left panel) and 
12+log(0/H) (right). The figure shows only those objects with 3cr detections or better in 
the [Oiv] line. 
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Fig. 7.— [Ne in]/ [Neil]/ and [Siv]/[Siii] (18.7/xm) plotted against the [0 iv] / [Si n] flux 



ratio s. As in Fig. H] BCDs are shown as filled circles, and the SINGS galaxies (IDale et al. 



20091 1 as open squares (Hll region-nuclei) and open triangles (AGN). The figure shows only 



those objects with 3a detections or better in the [O iv] line. 
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Fig. 8. — (z / L i/ ) 2 4/((^L v )7i+ (iaLv)i 60 (upper panel) and 12+log(0/H) (lower) plotted against 
the MIPS ratio ( yL u ) 7 -\ / (vL,,)-\ m . As in Fig. HJ BCDs are shown as filled circles, and the 
SINGS galaxies (IDale et al.l 120091 ) as open squares (Hll region-nuclei) and open triangles 
(AGN). Both panels show correlations significant within the BCD sample at ^ 99% confi- 
dence level, as described in the text. 



1 10 

<i/L„> 71 /<i/L„> 160 



Fig. 9. — [Ne ill] /[Nell]/ (lower panel) and [Olv]/[Sill] (upper) plotted against MIPS ratio 
(vL v ) 7 i I (tAL,;) 1fin . As in Fig. HJ BCDs are shown as filled circles, and the SINGS galaxies 
( Dale et al. 2009h as open squares (Hll region-nuclei) and open triangles (AGN). In the upper 
panel, the 3<r [Si n] upper limit (with [O iv] detection) is shown as an open circle. 



-54- 




io 2 io 3 

n e ([Sn]0.6717^m/[SII]0.6731^m) 



Fig. 10. — Electron densities n e derived from IR [S ill] line ratios vs. those inferred from 
optical [S ill] line ratios. The dotted line denotes equality. 
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Fig. 11. — [Si n] emission normalized by and plotted against TIR. The blue dotted horizontal 
line indicates the BCD mean, without any MIPS ULs and without UM311 (Ltir~ 3 x 
Q ), since its TIR is overestimated because of crowding (see tex t). The red dashed line 
corresponds to the mean for the SINGS sample (IRoussel et al.l 120071 ). 
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Fig. 12. — PAH strengths relative to total PAH luminosity vs. the total PAH power normal- 
ized to TIR. The four panels show the the following dust features (DF): 6.2 /zm (bottom left), 
7.7 /im (bottom right), 8.6 /im (top left), 11.3 (top right) /im features. BCDs are plotted as 
filled (blue) circles, with TIR lower limits shown as open circles; the TIR lower limits do not 
affect the vertical axis, only the placement along the horizontal one. Open (gr een) squares 



corre spond to SINGS HII nuclei, and open (red) triangles to SI NGS AGN (jSmith et al. 



20071 ). The horizontal dashed lines give the SINGS sample medians (jSmith et al.ll2007l ). and 
the dotted ones the means for the BCD sample, taking into account all objects with 7.7 /xm 
detections. The vertical dotted line corresponds to the BCD mean PAH power normalized 
to TIR. 
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Fig. 13. — Median Drude profile widths (FWHM) and central wavelengths for seven aromatic 
features detected ( £ 3a) in our sample. The error bars are the standard deviations of the 
measurements, and the numbers in parentheses are the numb ers of BCDs with significant 
detections. Open squares correspond to SINGS sample means (ISmith et al.l 120071 ) . 
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Fig. 14. — The integrated flux in the PAH features at 6.2 /im, 11.2 /im, and 12.7 /xm, plotted 
in two sets of ratios: 6.2/11.2 vs. 12.7/11.2. Only those BCDs with at least one ratio with 
11.2 /im are shown: filled (blue) circles correspond to those BCDs with 12+log (Q/H)>8.1 
open o nes to those objects with lower metallicities. The other data are taken from lHony et al 
(120011 ) and correspond to planetary nebulae (PNe, open red triangles), reflection nebulae 
(RNe, open red squares), intermediate- mass star-forming regions (six-sided open green stars), 
and Hll regions (open green hexagons). The SINGS aver age is shown as a (black) x, with 
error bars reporting standard deviations over the sample (ISmith et al.l 120071 ). 
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Fig. 15. — The integrated flux in all the PAH features detected at > 3a, relative to the total 
IR flux, plotted against the the total IR flux as described in the text. BCDs are plotted 
as filled blue circles. The non-AGN SINGS galaxies are plotted as open green squares, and 
the AGN SINGS galaxies (Seyferts and LINERs) as open red triangles. The blue dotted 
line shows the median S(PAH)/TIR for our sample of 0.54%; the red d ashed line shows the 
mean S(PAH)/TIR for the SINGS sample of -10% jSmith et al.ll20Q7h . 
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Fig. 16. — The integrated flux in all the PAH features detected at > 3a, relative to the total 
IR flux, plotted against the the oxygen abundance (left panel) and the [Ne ill] /[Nell] flux 
ratio (right). BCDs are shown as filled blue circles, the non-AGN SINGS galaxies as open 
green squares, and the AGN SINGS galaxies (Seyferts and LINERs) as open red triangles. 
The horizontal dotted line shows the BCD sample mean S(PAH)/TIR of 0.54%, and the 
dashed line the SINGS mean of 10%. The right panel shows only those objects with 3cr 
detections or better in each line; the best-fit linear regression is shown as a dotted line. 
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Fig. 17. — uL v (7. 7/jlWl) /TIR plotted again Ltir for the local starbursts in lHouck et al.l (120071 ) 
(shown as open red squares) and the BCDs (filled blue circles). The dotted horizontal line 
segments correspond to the sample means: 0.044 for the BCDs, and 0.167 for the starbursts. 



-61 - 



o 
u 



a. 

CO 



o 
o 



3. 




0.1 1.0 10.0 

F([NeIIl] 15.6/Lim)/F([NeIl] 12.8/zm) 



Fig. 18. — The ratios of the 7.7/11.3 /im bands plotted against [Ne ill] /[Nell]. BCDs with 
£ 3cr detections are shown as filled blue circles, SINGS AGN as open red triangles, and 
SINGS HII nuclei as open green squa res. The mean 7.7 /11.3 /xm ratio for the SINGS sample 
is shown as a horizontal dotted line (ISmith et all 120071 ) . 
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Fig. 19. — [Oiv][Neii] (or 1.7x33.5/xm [S in]) vs. the strength of the 7.7/xm PAH fe ature. 
The data for starburst galaxies, ULIRGs, and AGN are taken from lGenzel et al.l (119981 ). and 
are marked by green open stars, green open squares, and red open triangles, respectively. 
Our BCDs as filled blue circles, with upper limits for the PAH features shown by left-pointing 
arrows; only the (8) BCDs with £ 3a [O iv] dete c tions are shown. The "mixing curve" shown 
by a dashed line is made following iGenzel et al.l (119981 ); various fractions of total luminosity 
in an AGN and a starburst are combined, assuming a pure AGN has [O lv]/[Ne n]/ ~1, PAH 
strength ~0.04, and a pure starburst [Oiv][Neii]/ ~0.002, PAH strength ~3.6. The 50% 
mark is shown as a line segment along the mixing curve. 
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Fig. 20. — Diagnostic diagram as in iPeeters et al.l ( ]2004al ). adapted from lLaurent et al. 
(120001 ). The templates for the "mixing" curves are empirical, defined to be the positions 
of M17 (Hn region), NGC7023 (an exposed PDR), and a pu re AGN (assumed to have 
negligible 6.2 /im flux in the PAH band). Data are taken from IPeeters et al.l (l2004al ). and 
include Seyfert Is (filled red triangles), Seyfert 2s (open red triangles), ULIRGs (open green 
squares), starburst galaxies (open green stars), and Galactic star- forming regions (green x). 
Again, BCDs are shown as filled blue circles, with upper limits for the 6.2 /xm PAH feature 
shown as left-pointing arrows. 
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Fig. 21. — Sum of S(0) to S(2) H 2 luminosity normalized by and plotted against TIR. 
All detections with £ 2cr are plotted; a x denotes the 3a detections. Different symbols 
distinguish which H 2 lines are considered in the sum: Mrk996 with S(0) and S(l), and 
CGCG005 — 027 with S(l) and S(2), are shown by open circles; the remaining BCDs with 
only S(l) or S(2) by open squares. The red hatched area shows the range for SINGS galaxies 
( iRoussel et al.l 120071 ). and the blue dotted horizontal line indicates the BCD mean, without 
UM311 (L TIR ~ 3 x 1O 9 L ) because of the likely overestimate of its TIR because of source 
crowding (see text). 
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Fig. 22. — Excitation diagrams for the four BCDs with multiple significant H 2 detections 
and at least one significant ( £ 2.5a for SBS 1152+579) detection in a low-order transition 
(see text). 3a upper limits are shown with vertical arrows. The dotted lines correspond 
to best-fit estimates of the temperature of the line-emitting gas (see text). The parabolic 
fit for CGCG 005 — 027 is an experiment to show the continuity of the data, and thus the 
continuous range of gas temperatures. 
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